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The  purpose  of  this  study  is  to  determine  how  wind  shear 
affects  the  lateral  spread  of  fallout  patterns.  Several 
mathematical  expressions  are  currently  in  use  to  predict  this 
growth.  Determining  which  one  is  the  most  accurate  is  the  focus 
and  goal  of  this  study.  It  may  even  be  that  none  of  these  are 
correct.  In  that  case,  a  new  expression  will  be  sought  and,  if 
possible,  validated. 

Fallout  patterns  clearly  extend  in  both  the  downwind  and 
crosswind  directions.  In  the  crosswind  direction,  perpendicular 
to  the  downwind  direction,  the  fallout  is  usually  considered  to 
be  distributed  in  a  gaussian,  or  near-gauss ian,  manner.  This 
gaussian  distribution  can  be  described  by  its  characteristic 
width,  o  .  Thus,  a  is  a  good  predictor  of  the  total 
lateral  distribution  of  the  fallout  pattern.  The  models  used 
throughout  this  study  are  oriented  on  determining  this  charac¬ 
teristic  o  . 
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Abstract 

This  study  refines  the  method  of  determining  the  lateral 
spread  of  fallout  from  nuclear  explosions.  Lateral  growth  of  a 
fallout  "pattern"  under  the  influence  of  torroidal  growth,  wind 
shear,  and  diffusive  growth  is  examined  to  determine  how  best  to 
calculate  the  extent  of  a  fallout  pattern  as  it  grows  with  time. 

Two  models  of  the  initial  stabilized  nuclear  cloud  are 
investigated  to  determine  how  best  to  model  the  initial 
conditions  which  define  the  starting  positions  of  the  fallout 
particles.  A  hybrid  model,  incorporating  portions  of  both,  is 
developed.  The  grounded  locations  of  the  fallout  are  calculated 
after  examining  the  influence  of  horizontal  winds  on  their 
trajectories.  From  this  examination,  with  the  assistance  of  a 
new  mathematical  derivation  of  o s,  an  improved  formula  for  the 
prediction  of  lateral  spread  of  fallout  is  tested  and  confirmed. 

Certain  fairly  common  assumptions  about  the  atmosphere 
and  the  use  of  the  DELFIC  100  e qu a  1 -ac t i v  i  t y  groups  allow 
simplification  of  the  treatment  of  the  geophysical  properties  to 
the  extent  that  the  problem  can  Ke  analyzed  on  a  microcomputer 
using  BASIC.  Accuracy  is  not  sacrificed  for  simplicity,  yet  the 
methods  used  can  provide  quicker,  more  accurate  results  than  the 
large  codes  run  on  mainframe  computers. 

Although  only  a  one  megaton  explosion  is  examined  in 
detail,  the  methods  are  applicable  to  the  full  range  of 
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strategic  nuclear  weapons. 


DETERMINATION  OF  LATERAL  SPREAD 


OF  NUCLEAR  FALLOUT  PATTERNS 


I.  Introduction 


More  than  ever  before,  accurate  and  timely  information 
has  become  the  cornerstone  of  sound  tactical  and  strategic 
decisions  by  military  leaders.  In  this  age  of  sophisticated 
nuclear  weapons  and  near  "real-time"  intelligence,  brought  about 
by  micro-miniaturization  and  advanced  electronics,  it  is 
surprising  to  find  that  little  has  been  done  to  provide  planners 
with  methods  of  determining  nuclear  weapons  effects  that  are 
al3o  advanced,  quick,  and  relatively  easy  to  use.  In  fact,  much 
of  today's  planning  is  based  on  research  that  was  stopped  in  the 
1960's.  Although  atmospheric  testing  of  nuclear  weapons  was 
stopped  as  a  result  of  the  Nuclear  Test  Ban  Treaty,  there  remain 
areas  of  ambiguity  and  uncertainty  in  the  analysis  of  the  data 
obtained  during  those  tests. 

Advances  in  peripheral  areas,  such  as  atmospherics,  can 
be  applied  to  the  original  data  from  atmospheric  tests  to  pro¬ 
vide  greater  understanding  of  the  processes  involved.  Increased 
understanding  could  lead  to  improvements  in  accuracy  and  speed. 


Although  some  nuclear  effects  may  never  be  fully  understood, 


others  can  already  be  more  fully  explained  with  the  current 
state  of  knowledge. 

Hopefully,  this  study  will  further  understanding  in  one 
area  of  nuclear  weapon  effects.  That  area  is,  specifically,  the 
growth  of  the  lateral  spread  of  the  fallout  pattern  with  respect 
to  time  due  to  crosswind  shear.  From  an  analysis  of  the  forces 
which  create  nuclear  clouds,  the  transportation  of  those  clouds 
through  the  atmosphere,  and  the  resulting  pattern  of  fallout 
created  on  the  ground,  an  expression  relating  the  time  of 
arrival  of  the  fallout  and  the  characteristic  width  of  the 
fallout  pattern  (defined  as  "sigma",  o  ,  from  the  standard 
definition  of  a  Gaussian  distribution)  can  be  developed  for  a 
variety  of  initial  conditions.  The  expressions  developed  in  this 
manner  will  be  compared  against  the  optimum  "full-physics"  model 
of  fallout  pattern  formation  to  determine  which  expression  most 
closely  replicates  the  actual  process. 


II.  Problem  Statement 


There  are  currently  two  primary  equations  which  compete 
for  use  in  determining  the  lateral  growth  of  fallout  patterns. 
These  equations,  which  are  explained  in  the  next  chapter,  are 
both  empirical  or  semi-empirical  and  are  the  result  of  different 
methods  of  analysis,  yet  they  are  supposed  to  predict  the  same 
degree  of  lateral  growth  in  a  particular  situation  when  they  are 
used.  This  research  effort  focuses  on  determining  which  of  these 
aquations  more  accurately  predicts  the  lateral  growth.  There 
exists  the  possibility  that  neither  one  of  them  predicts  the 
fallout  pattern  as  well  as  a  new  and  different  equation  which  is 
presented  here. 

In  short,  the  purpose  of  this  research  is  to  determine 
what  equation  best  predicts  the  lateral  growth  of  a  fallout 


pattern  with  respect  to  time. 


III.  Background 


Unique  Terms 

Fallout  prediction  uses  some  terms  that  can  be  confusing 
when  their  specific  definition  within  the  field  is  not  made 
clear.  As  a  starting  point,  five  definitions  are  included  here. 

Downwind.  It  is  convenient  to  use  a  standard  cartesian 
coordinate  system  as  an  initial  reference  frame.  As  can  be  seen 
in  Figure  1,  the  X  axis  is  co-labeled  as  the  "Downwind"  axis. 
Because  the  process  of  fallout  formation  occurs  in  the 
atmosphere,  it  is  critical  to  define  the  orientation  of  this 
coordinate  system  with  respect  to  that  same  atmosphere.  For  the 
purposes  of  this  research,  the  following  convention  is  chosen 
for  determining  the  orientation  of  the  downwind  axis.  The  wind 
direction  at  the  altitude  of  the  center  of  the  stabilized 
nuclear  cloud  (this  term  will  be  further  defined  in  Chapter  V) 
is  identical  to  the  direction  of  the  X  axis.  The  downwind 
direction  is,  then,  the  direction  in  which  a  nuclear  cloud  would 
drift  if  the  wind  were  blowing  from  the  same  direction  at  all 
altitudes . 

Cro3swind  .  Just  as  the  Y  axis  is  translated  ninety 
degrees  from  the  X  axis  in  a  cartesian  coordinate  system,  the 
Cros3wind  axis  represents  an  axis  that  has  been  rotated  ninety 
degrees  from  the  downwind  axis.  Using  this  definition, the  cross- 
wind  term  refers  to  any  component  of  the  wind  or  any  other  force 
which  acts  perpendicular  to  the  downwind  axis.  (See  Figure  1) 


Shear.  "Shear"  is  defined  as  the  magnitude  of  the  change 
in  both  speed  and  direction  between  any  two  wind  vectors. 
Naturally,  all  wind  vectors  occuring  in  the  atmosphere  will  have 
three  components  corresponding  to  the  three  axes  of  the 
reference  frame.  The  specific  formulation  of  the  shear  forces  as 
they  pertain  to  this  study  are  further  defined  in  Chapter  VI. 

Sigma.  The  standard  definition  of  a  gaussian,  or  Normal 
probability  function  (1:927)  is, 

1  -1/2[(x-u)/o]2 


where  the  probability  of  an  occurrence,  P(x),  is  defined  as  a 
function  of  the  mean,  y  ,  and  the  standard  deviation,  o  .  "Sigma", 
or  o  ,  is  alternatively  defined  as  the  standard  deviation,  the 
dispersion,  or  even  as  the  characteristic  width  of  the 
distribution.  The  fallout  pattern  on  the  ground  will  be  modeled 
with  the  gaussian  distribution  in  the  crosswind  direction.  Why 
this  is  the  case  will  be  shown  later.  Therefore,  oy  is  the 
characteristic  width  of  the  fallout  pattern  in  the  crosswind 
direction . 

Di3ktosser.  A  "Disktosser"  is  the  common  term  for  a  type 
of  numerical  fallout  code  which  uses  a  discrete  mesh  of  mono¬ 
sized  particle  disks  to  model  the  nuclear  cloud.  The  disktosser 
"tosses"  each  of  the  disks  that  characterize  the  initial  cloud 
into  a  given  wind  field  and  computes  a  grounded  location  for 
each  disk  based  on  the  wind  field  and  the  parameters  of  the 


Unfortunately,  there  were  simply  too  many  external  variables 


which  were  not  recorded  during  the  various  tests.  An  example  of 
the  difficulty  of  such  an  approach  might  best  be  illustrated 
by  a  quotation  from  The  Effects  of  Nuclear  Weapons  by  S. 
Glasstone  and  P.J.  Dolan  where  they  explain  a  fallout  pattern 
from  a  megaton-size  explosion. 


"However,  other  patterns  are  possible;  one,  for 
example,  ascribes  the  large  radiation  doses  on  the 
northern  islands  of  Rongelap  Atoll  to  a  hot  spot  and 
brings  the  3000-rad  contour  line  in  much  closer  to 
Bikini  Atoll.  Because  of  the  absence  of  observations 
from  large  areas  of  ocean,  the  choice  of  the  fallout 
pattern  is  largely  a  matter  of  guesswork.  (2:436) 


For  the  most  part,  nuclear  tests  conducted  over 
land,  although  not  in  the  megaton  range,  are  also  a  product  of 
the  same  kind  of  analysis.  In  short,  there  are  too  many 
variables  to  conduct  a  thorough  a  posteriori  analysis  of  past 
weapons  tests  at  the  moment. 


The  second  method,  which  relies  less  on  published  data 
and  more  on  the  physics  of  the  actual  process  of  transporting 
the  particles  through  the  atmosphere,  becomes  the  method  of 
choice  after  reviewing  the  available  published  data  on 
atmospheric  shots.  This  method  involves  modeling  the  individual 
phenomena  as  closely  as  possible  with  computer  models.  Each 
model  is  designed  to  duplicate,  as  closely  as  possible,  the 
actual  processes.  (  Such  codes  might  be  called  "full-physics" 
codes.)  These  might  be  divided  into  three  main  phases  of  fallout 
pattern  formation.  These  three  phases  are  1)  cloud  formation,  2) 


particle  fall,  and  3)  particle  accumulation  on  the  ground. 
Combining  all  three  phases  produces  the  location,  on  the  ground, 
of  a  particular  crosswind  portion  of  the  complete  grounded 
pattern.  Adding  all  these  portions,  or  slices,  results  in  a 
pattern  that  can  be  analyzed  to  determine  how  the  pattern  grew 
with  time.  Then,  the  pattern  of  growth  predicted  by  the  model 
will  be  compared  with  the  two  expressions  to  be  shown  next. 

The  major  disadvantage  of  the  modeling  method  is  that  it 
does  not  lend  itself  to  comparison  with  an  actual  explosion. 
However,  its  ability  for  rigid  comparison  between  the  parameters 
and  analysis  of  the  physics  involved  allow  for  more  precise 
isolation  of  the  wind  shear  term.  This  is  a  clear  advantage. 

The  Two  Primary  Equations 

WSEG- 1 0.  One  of  the  pioneering  efforts  to  quantify  and 
simplify  what  was  then  known  about  nuclear  fallout  is  WSEG 
Research  Memorandum  No. 10  -  A_n  Analytic  Model  of  Close-In 
Deposition  of  Fal lout  for  Use  in  Operational -Type  Studies.  This 
report,  written  in  1959,  has  been  modified  several  times  in  the 
intervening  years  in  an  attempt  to  more  closely  match  the 
experimental  data  obtained  from  the  atmospheric  tests  (3:3-^). 
However,  the  method  of  calculating  the  lateral  growth  of  the 
fallout  pattern  has  not  been  significantly  modified  from  the 
original  method  as  shown  in  WSEG-10.  The  only  modification  to 
the  original  formulation  is  a  refinement  of  the  low-yield  (in 
the  kiloton  range)  deposition  rate.  This  refinement  produces  an 
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enhanced  fallout  pattern  in  the  vicinity  of  ground  zero  for 
early  times.  The  change  is  achieved  by  refining  the  input 


parameters  (4:2-3),  though,  not  by  altering  the  basic  method  of 
calculation.  Therefore,  the  underlying  formulation  of  o  as  used 
by  WSEG-10  can  be  extracted  from  their  original  equation  which 
defines  the  total  crosswind  spread.  The  total  crosswind  spread 
is  caused  by  three  main  contributors,  which,  depending  on  the 
time  at  which  the  fallout  pattern  is  analyzed,  have  different 
degrees  of  effect.  The  three  contributors  are  cloud  rise 
spreading,  torroidal  spreading  after  cloud  stabilization,  and 
crosswind  shear  effects.  The  complete  formulation  of  crosswind 
spread  (5:14)  is, 

Oy2  =  o02  +  [(8xoo2)/L]  +  [SyOh^/W)  ]2  (km)  (1) 

where 


0y  is  the  total  crosswind  spread  (km) 

Og  is  the  width  of  the  initial  stabilized  cloud  center  (km) 
x  is  the  downwind  distance  at  which  the  spread  is  desired 
L  is  an  empirical  constant  determined  by  yield 
Sy  is  the  effective  shear  (km/km-hr) 

oh  is  the  width  of  the  vertical  distribution  of  the 
stabilized  cloud  (km) 

W  is  the  effective  wind  speed  in  the  downwind  direction. 
Since  x=Wta  where  ta  is  the  arrival  time,  the  WSEG  form 


may  be  written  as, 


(km) 


(la) 


oy2=o02+8ta°02/Tc+^syta3h)2 

where  TC=L/W  is  also  empirical  and  the  other  variables  are 
defined  above. 

With  the  aid  of  an  assumption,  equation  (1)  can  be 
reduced  as  follows.  The  second  term,  (8xoq^/L),  representsthe 
growth  due  to  torroidal  circulation.  As  a  result,  of  obser¬ 
vations  of  actual  clouds,  WSEG  suggests  that  this  growth 
continues  for  the  first  three  hours  of  cloud  formation  and, 
thereafter,  ceases  to  be  significant.  Thus  for  ta  greater  than  3 
hours  the  second  term  also  becomes  a  constant.  And, 

Oy2  ( t )  =  og^+1*oo^/Tc  +  (SyOi1ta  )  2  =  (  o^*  )  2  +  ( Sy  Oftta  )  2  (km)  (2) 

where  03*  is  the  cloud  width  at  3  hours  less  the  shear  growth 
effects. 

Without  derivation,  WSEG  (Weapon  Systems  Evaluation 
Group)  determined  that  the  shear  contribution  to  lateral  growth 
was  the  final  term  (Sy0^ta).  Designating  this  term  by  os,  and 
solving  for  os  produces  an  equation  that  defines  the  spreading 
due  to  shear  according  to  WSEG  as, 

°s=V°y2"  °3*2  (km)  (3) 

The  authors  do  not  include  a  derivation  of  equation  (1) 
anywhere  in  WSEG-10  (5:11).  However,  its  derivative  (equation 
(3))  will  still  serve  as  a  base  line  against  which  other 
expressions  for  the  growth  of  0  will  be  compared. 


10 


Norment.  The  second  major  formulation  of  the  growth  of  o 
which  has  achieved  a  degree  of  acceptance  is  put  forward  by  H. 
G.  Norment  in  "DNAF-1:  An  Analytical  Fallout  Prediction  Model 
and  Code.”  He  presents  a  form  similar  to  WSEG-10's  formulation 
with  replaced  by  the  total  cloud  vertical  thickness  and  the 
addition  of  a  divisor  of  10.  Norment's  formula  is  (6:41), 

°s=[sy (zT-ze)ta]/10  (km)  (4) 

where 


is  the  stabilized  cloud  top, 
zq  is  the  stabilized  cloud  bottom, 
ta  is  the  time  of  arrival  of  the  fallout. 

Norment  states  that  this  equation  is  "somewhat 
arbitrary,"  but  "chosen  by  numerical  experimentation  to  give 
good  comparisons...."  (6:44)  Equation  (4)  is  taken  as  the  second 
empirical  prediction  of  the  growth  of  the  lateral  spread  of  the 
fallout  pattern  against  which  to  compare  the  results  of  this 
study. 

Approach 

The  development  of  the  assumptions  necessary  to  establish 
a  base  line  are  covered  in  Chapter  IV. 

Chapter  V  discusses  nuclear  cloud  formation  as  presented 


by  WSEG-10  and  Hopkins.  It  also  compares  the  two  models  and 
finds  the  most  appropriate  one  to  use  for  this  study. 


The  method  of  modeling  the  fall  mechanics  as  they  apply 
to  micron  sized  particles  is  addressed  in  Chapter  VI.  This 
chapter  also  includes  a  description  of  the  "full-physics" 
computer  code  used  as  a  benchmark. 

Chapter  VII  shows  how  the  fallout  deposited  on  theground 
is  treated  to  determine  the  characteristic  width. 

Chapter  VIII  compares  the  Norment  and  WSEG  equations  to 
the  spread  predicted  by  the  model  and  presents  a  new  derivation 
by  Bridgman  which  is  also  compared. 

Conclusions  and  recommendations  are  presented  in  Chapters 


IX  and  X 


IV.  Assumptions 


In  order  to  establish  a  starting  point  for  the  study, 
several  assumptions  are  made.  This  chapter  details  those 
assumptions  and  justifies  them. 

1 )  Local  Fallout 

Local  fallout  is  hereby  specifically  defined  as  that 
fallout  which  falls  to  the  earth  within  the  first  twenty-four 
hours  following  detonation  time.  This  definition  is  somewhat 
more  restrictive  than  some.  However,  as  will  be  shown,  one  of 
the  major  areas  of  controversy  is  how  the  sub-micron  size 
particles  interact  with  the  atmosphere.  Precisely  whether  they 
fall  by  gravity  or  are  carried  by  the  turbulent  atmosphere,  sub¬ 
micron  and  even  few  micron  particles  are  not  grounded  in  24 
hours  (7:Ch  2,16).  This  controversy  can  be  minimized  by  limiting 
the  calculations  to  the  first  twenty-four  hours. 

It  is  further  assumed  that  all  of  the  fallout  that  is 
present  in  the  initial  stabilized  cloud  can  be  described  by  the 
100  equal-activity  groups  determined  by  Conners  as  derived  from 
the  DELFIC  (Defense  Land  Fallout  Information  Code)  computer  code 
(8:14).  Furthermore,  any  of  the  particles  thus  described  which 
are  determined  to  fall  to  the  ground  within  twenty-four  hours 
will  not  have  been  influenced  by  any  buoyant  forces  inherent  in 
the  atmosphere.  In  summary,  the  particles  of  fallout  are  modeled 


as  spheres,  with  equal  mass  radii  determined  by  the  DELFIC 
default  spectrum,  subject  to  aerodynamic  drag  as  they  fall 
through  the  atmosphere. 

2 )  Stratified  Atmosphere 

The  earth's  atmosphere  is  a  complex  region  where  many 
forces  interact  concurrently,  from  the  earth's  surface  through 
the  ozonosphere,  to  create  air  movement.  However,  even  the  most 
powerful  thermo-nuclear  weapons  (300  megatons)  are  not  thought 
to  be  capable  of  creating  clouds  that  will  stabilize  above  the 
stratosphere.  The  troposphere  and  the  stratosphere,  then,  are 
the  regions  of  the  atmosphere  which  are  of  most  concern  for  this 
study.  Within  these  regions,  the  temperature,  density,  pressure, 
and  dynamic  viscosity  are  assumed  to  be  as  defined  in  U.S. 
Standard  Atmosphere  (9). 

The  winds  in  the  troposphere,  being  easier  to  measure 
and  of  more  immediate  concern  to  many  researchers,  are  fairly 
well  documented.  Changes  in  wind  speed  and  direction  occur 
gradually  from  point  to  point.  The  magnitude  of  the  change  in 
wind  speed  and  direction  with  altitude  is  defined  as  horizontal 
shear.  Numerical  methods  of  analysis  must  represent  this  smooth 
change  by  discontinuous  approximations.  Available  data,  as  is 
presented  in  U.S.  Standard  Atmosphere,  is  usually  in  discrete 
elements.  U.S. Standard  Atmosphere  uses  50  meters  below  11.1  Km 
and  100  meters  for  the  region  between  11.1  Km  and  20.0  Km.  But, 


when  the  magnitudes  of  shear  winds  can  be  found,  they  are 


usually  reported  in  1000  meter  increments.  As  a  starting  point, 
it  is  assumed  that  a  vertical  increment  of  200  meters  gives  an 
adequate  representation  of  shear  for  the  purpose  of  this 
analysis.  This  200  meter  increment  is,  admittedly,  a  compromise 
between  the  degree  of  accuracy  available  by  using  the  same  fifty 
meters  as  U.S.  Standard  Atmosphere  and  the  speed  of  computation 
which  comes  from  using  an  increment  oflOOO  meters. 

The  stratosphere,  by  contrast,  is  considerably  less  well 
documented.  Even  though  wind  currents  appear  to  be  present  at 
all  altitudes  within  the  troposphere  and  the  stratosphere,  it  is 
often  customary  to  consider  the  stratosphere  to  be  composed  of 
very  large  air  masses  moving  essentially  geostrophical ly.  These 
geostrophic  winds  would  negate  shear  except  at  the  boundaries  of 
the  air  masses.  Thus, some  analyses  involving  the  stratospheric 
winds  disregard  shear  forces  completely  (10:81).  However,  since 
some  air  movement  in  the  stratosphere  can  be  documented  and  the 
precise  lower  boundary  of  the  stratosphere  changes  from  day  to 
day,  this  analysis  will  assume  that  there  is  some  degree  of 
shear  at  all  altitudes  of  interest.  The  magnitude  and  direction 
of  shear  at  all  altitudes  is  adressed  in  greater  detail  in 
Chapter  VI. 

3 )  Uniform  Flat  Earth 

Inclusion  of  terrain  features  in  these  kinds  of  models 
causes  each  model  to  predict  results  which  are  uniquely  defined 
by  the  type  and  size  of  the  included  terrain.  This  makes 


comparison  between  different  studies  very  difficult.  Most  such 
studies  assume  a  uniform  flat  earth. 

Additionally,  treating  the  earth's  surface  as  a  smooth 
surface  allows  the  model  to  continue  calculating  wind  shear 
effects  linearly  until  downward  movement  of  the  fallout  has 
stopped.  There  are  two  alternatives  to  this  method.  The  first 
alternative  is  to  incorporate  the  boundary  layer  turbulence 
effects  for  the  viscous  sublayer  (10:155).  Such  a  treatment  is 
difficult  and  is  restrictive  to  a  specific  location.  The  end 
result  is  a  shortened  fall  distance  and  a  smaller  total  fallout 
pattern.  The  other  alternative  is  to  simply  ignore  continued 
shear  effects  in  the  region  close  to  the  surface.  This  results 
in  simplifying  the  calculation  but  shortens  the  actual  height  of 
fall.  The  result  is  the  same  as  if  the  first  alternative  is 
used . 

Therefore,  the  models  used  in  this  study  assume  a  uniform 
flat  earth  with  a  constant  gravitational  field.  It  is  recognized 
that  this  assumption  causes  the  resultant  fallout  pattern  to  be 
larger  than  if  surface  roughness  is  included.  This  increase  in 
the  pattern  is  addressed,  and  accounted  for,  during  phase  three 
(Chapter  VII). 

4 )  Constant  Activity 


Actual  fallout  patterns  incorporate  radioactive  decay  of 
the  fallout  particles  during  their  fall.  This  is  an  important 
feature  in  calculating  radiation  doses  on  the  ground.  For  the 


purpose  of  determining  where  the  activity  is,  though,  the 
actual  amplitude  of  that  activity  is  unimportant.  By  working 
with  unit  time  reference  activity,  (the  activity  at  1  hour  after 
burst  =  530E1C)6  gamma  curies/Kt),the  pattern  from  point  to  point 
can  be  compared  directly  (2:453). 

5 )  Shear  Dominates  Only  After  Three  Hours 

Torroidal  circulation  within  the  expanding  cloud  is 
believed  to  be  the  dominant  effect  that  continues  to  expand  the 
cloud  from  stabilization  time  until  approximately  three  hours 
after  burst  time.  This  process,  although  not  well  understood,  is 
treated  by  observations  of  the  actual  clouds.  No  attempt  is  made 
to  calculate  the  effects  of  torroidal  growth  in  the  model  which, 
thus,  does  not  examine  the  period  from  detonation  to  three 
hours.  It  is  assumed  that  torroidal  circulation  has  a  negligible 
effect  on  cloud  growth  from  three  hours  on. 


The  assumptions  provide  the  initial  framework  within 
which  to  begin  the  analysis  of  cloud  growth.  This  chapter 
describes  the  stabilized  cloud.  Because  subsequent  full-physics 
codes  use  the  stabilized  cloud  as  the  initial  conditions  for 


computing  the  fall  of  the  individual  particles,  it  is  necessary 
to  model  the  initial  cloud  as  completely  as  possible. 

Approximately^ i ve  to  seven  minutes  after  a  surface  burst 
of  a  one  megaton  nuclear  device,  the  majority  of  the  upward 
movement  of  the  particulate  matter  caught  up  by  the  fireball  and 
afterwinds  is  considered  to  have  slowed  sufficiently  that  the 
earth's  gravity  begins  to  pull  the  particles  back  down. 


Torroidal  circulation  within  the  vortex  of  the  cloud  continues 


to  affect  the  location  of  the  particles  within  the  cloud  to  some 
extent.  This  circulation  causes  continued  expansion  in  both  the 
vertical  and  horizontal  directions  until  approximately  three 
hours  after  burst.  However,  the  dominant  effect,  which  causes 
the  upward  movement  to  stop,  is  the  cooling  of  the  material  in 
the  cloud  to  such  a  degree  that  it  essentially  reaches  an 
equilibrium  state  with  the  surrounding  air.  This  rise  is  heavily 
dependent  on  the  altitude  of  the  tropopause  and  the  heat  energy 
remaining  in  the  cloud  at  altitude  (2:31-32).  The  upward 
movement  stops  while  the  horizontal  torroidal  expansion 
continues.  The  point  in  time  when  this  vertical  equilibrium 
state  is  reached  is  called  the  stabilization  time. 

Compared  to  the  hours  which  are  required  for  much  of  the 


activity  to  settle  to  the  ground,  the  initial  six  to  seven 
minutes  of  cloud  rise  time  seem  to  be  negligible.  Unless  very 
strong  winds  prevail  throughout  the  short  rise  time  of  the 
cloud,  the  cloud  rises  relatively  straight  upward  along  a 
vertical  axis.  This  study  does  not  look  at  wind  conditions  that 
are  strong  enough  to  cause  significant  movement  off  this  axis. 
This  allows  the  initial  stabilized  cloud  to  be  modeled,  as  shown 
in  Figure  2,  with  the  center  axis  of  the  cloud  directly  over 
the  burst  point  (also  called  Ground  Zero).  The  shape  of  the 
cloud  is  modeled  as  an  upright  circular  cylinder  that  resembles 
a  tin  can.  For  the  moment,  the  actual  dimensions  of  the  cloud 
remain  unspecified. 


Figure  2.  The  Initial  Cloud  Model 


WSEG's  Stabilized  Cloud  Dimensions 


During  the  writing  of  WSEG-10,  the  authors  noted  that  by 
plotting  the  activity  present  in  the  cloud  versus  its  altitude, 
as  determined  from  NCTS  rocket  data,  the  spatial  distribution  of 
activity  could  be  adequately  represented  by  three  orthogonal 
distributions  (5:17-19).  The  three  distributions  corresponded 
with  the  three  standard  cartesian  coordinate  axes.  By  examining 
each  of  the  three  distributions  separately,  they  found  that 
although  there  were  some  regions  of  poor  fit,  (particularly  for 
the  lateral  distribution),  that  there  was  good  enough  agreement 
to  consider  each  distribution  to  be  gaussian  (5:21). 

The  "effective"  cloud  radius  was  then  defined  as  "the 
initial  width  of  the  pattern"(5: 19)  and  then  by  comparison  with 
D.A.S.A.  data  for  the  "cloud  radius  normalized  to  fit  ground 
zero  dose  rate"  (5:24),  they  were  able  to  arrive  at  the 
following  approximation  for  determination  of  the  characteristic 
width  of  the  effective  initial  cloud  radius. 


lno0=0.70+ln(Y)/3-3.25/[4.00+(ln(Y)+5.4)2] 


where 


°0  is  in  statute  miles 
Y  is  the  yield  in  megatons 


Equation  (5)  has  yet  to  be  improved  upon  and  its 
improvement  is  not  contemplated  here. 

By  examination  of  the  distribution  of  activity  in  the 


vertical  direction,  G.Pugh  and  R.  Galiano  (the  WSEG-10  authors) 
obtained  a  curve  of  radioactive  cloud  center  height,  hG*  versus 
yield  that  could  be  fit  by  the  following  equation  (5:24), 

hc=44+6. 11n(Y)-.205[ln(Y)+2.42] |ln(Y)+2.42|  (6) 

where 

Y  is  the  yield  in  megatons 

hc  is  in  kilofeet. 

(To  avoid  confusion,  the  subscript  "0"  which  is  sometimes 
used  in  WSEG-10  to  indicate  the  "center"  is  here  replaced  by  the 
subscript  "c".) 

And  from  the  same  D.A.S.A.  data  they  concluded  that  the 
characteristic  vertical  width  at  hc  was, 

oz=0.18  hc  (7) 

where  oz  and  hc  are  in  kilofeet  (3:24). 

The  lateral  component  of  the  orthogonal  gaussian 
distribution  simply  spreads  the  total  activity,  defined  by  the 
amplitude  of  the  vertical  gaussian  at  a  given  altitude, 
horizontally  across  the  diameter  of  the  cloud  at  the  altitude 
for  which  the  amplitude  was  obtained.  As  with  the  other  two 
distributions,  torroidal  circulation  affects  the  size  of  the 
area  over  which  this  activity  fraction  is  spread.  WSEG-10  first 
U3es  the  assumption  that  torroidal  circulation  constricts  the 
initial  stabilized  cloud  width  to  about  half  of oq.  A  second 
assumption  that  the  torroidal  effect  diminishes  to  0,  at 


approximately  three  hours,  is  used  to  achieve  some  measure  of 
the  torroidal  growth  (5:13— However,  if  the  behavior  of  the 
fallout  pattern  near  the  origin  is  not  significant  and  if 
torroidal  growth  is  assumed  to  be  a  constant  contributor  to  the 
overall  width  of  the  pattern  after  three  hours,  then,  when  only 
times  after  three  hours  are  examined  for  a  determination  of 
growth  due  to  shear,  the  total  contribution  due  to  torroidal 
circulation  can  be  ignored. 

Figure  3,  below,  is  included  to  assist  visualization  of 
the  three  distributions  which  describe  the  initial  stabilized 
cloud  according  to  WSEG-10. 
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Figure  3.  The  WSEG-10  Model  Stabilized  Cloud 
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It  is  perhaps  interesting  to  look  at  what  equations  (5), 
(6),  and  (7)  produce  for  a  variety  of  weapon  yields.  Table 
I,  below,  lists  a  range  of  yields  and  the  values  WSEG-10 
predicts  for  those  yields. 


Table  I 

WSEG  Predicted  Stabilized  Cloud  Dimensions 

Yield 

MM 

°0 

(Mt) 

(Km) 

0. 17 

10.088 

1.816 

1.485 

1 

13.043 

2.348 

2.938 

5 

15.386 

2.769 

5.212 

9 

16.161 

2.909 

6.394 

20 

17.145 

3.086 

8.420 

100 

18.886 

3.399 

14.577 

Table  I  lists  a  100  megaton  yield.  As  of  this  writing,  it 
is  not  believed  that  such  a  large  warhead  is  in  the  inventory  of 
any  nation.  Using  the  temperature  change  as  a  guide  to 
determining  both  the  top  of  the  troposphere  and  the  top  of  the 
stratosphere  from  U.S.  Standard  Atmosphere,  one  can  see  that  the 
troposphere  ends  at  11.10  km  and  the  stratosphere  at  20.00  km 
(9:59,61).  Other  sources  allow  the  stratosphere  to  extend  to  as 
high  as  30  km  (11:194).  Howe ver , Tab  1 e  I  points  out  the  fact 
that  even  the  very  large  yield  warheads  are  still  predicted  to 
stabilize  within  a  conservatively  estimated  stratosphere. 
Similarly,  inspection  of  Table  I  shows  that,  conservatively,  a 
major  portion  of  the  radioactive  cloud  from  megaton  size  weapons 
will  be  above  the  troposphere.  It  is  because  of  the  uncertainty 
in  the  precise  limits  of  the  troposphere  and  stratosphere,  and 

3s 


the  near  certainty  in  stabilization  heights  as  predicted  by  the 


WSEG-10  equations  that  the  second  assumption  regarding  shear  in 
the  stratosphere  is  so  important  to  the  determination  of  os* 

Lastly,  Table  I  shows  the  relevance  of  using  a  one 
megaton  yield  for  this  study.  Namely,  it  represents  median  range 
values  for  yield,  stabilization  height,  and  cloud  radius. 

Hopkins  *  Stabilized  Cloud  Dimensions 

The  WSEG-10  model  relies  heavily  on  modeling  the  three 
spatial  distributions  as  smooth  gaussian  distributions.  The 
three  spatial  variables,  which  are  assumed  to  be  independent  of 
each  other,  can  be  considered  in  isolation  to  determine  whether 
or  not  they  can  be  improved  upon.  In  the  case  of  the  two  hori¬ 
zontal  distributions,  they  appear  satisfactory,  since  the  early 
effects  of  the  winds  are  discounted  during  cloud  formation,  and 
the  stratified  atmosphere  assumes  that  the  distributions  will 
encounter  the  changing,  non-homogeneous  medium  at  all  points 
along  the  distributions  at  the  same  moments  during  their  ascent. 
Nothing  should  cause  any  distortion,  then,  to  the  horizontal 
distributions  of  activity  during  cloud  rise.  The  vertical  dis¬ 
tribution  must,  however,  encounter  every  non-homogeneity  in  the 
atmosphere  at  succeeding  times  as  it  rises  and  spreads  out 
vertically.  This  gives  rise  to  the  speculation  that  a  more 
sophisticated  model  of  this  process  must  be  utilized. 

A.  T.  Hopkins  presented  in  his  thesis,  k  Two  Step  Method 
To  Treat  Variable  Winds  In  Fal lout  Smearing  Codes,  a  different 
method  for  modeling  the  vertical  distribution  of  activity  in  the 


stabilized  cloud  (12:  14-15).  Essentially,  the  method  considers 
each  of  the  particles  defined  by  the  DELFIC  100  equal-activity 
groups  individually  and  lofts  them  to  altitudes  that  are 
dependent  on  the  mass  of  each  particle.  This  method  determines  a 
separate  stabilization  height  for  each  activity  group.  A  unique 
oz  was  determined  by  Conners  for  the  distribution  of  each 
activity  group  in  the  vertical  direction  (8:19).  The  100 
separate  distributions  add  together  to  form  the  total  vertical 
distribution . 

By  running  the  DELFIC  code  with  yields  ranging  from  1  Kt 
to  15  Mt,  obtaining  a  fit  of  particle  size  versus  altitude  with 
a  linear  function  for  each  yield  as  an  intermediate  result,  and 
then  fitting  the  slopes  and  intercepts  of  the  intermediate 
results  to  polynomials  expressed  as  logarithms  of  the  initial 
yields,  Hopkins  (12:14-15)  determined  that 


where 


h0i  =  Im<-2rmlsm  <"> 

ho*  is  the  stabilized  altitude  of  the  ith  group, 

Im  is  the  intercept  of  the  radius  in  meters, 

rm*  is  the  mean  radius  of  the  it^  group  in  meters, 

Sm  is  the  slope  of  the  altitude  per  meter, 


and  both  Im  and  Sm  are  given  by 


Im=EXP(7.889+0.341n(Y>+  0 . 00 1 226  In ( Y )2 
-0.0052271n(Y)3+0.0004171n(Y)a) 


Sm=-EXP(1 .574-0. 011 971n(Y)+0.03636ln(Y)2 


-0. 0041  In (Y)3+0. 0001 9651n(Y)4) 


Table  II,  below,  provides  a  representative  sample  of  what 


Hopkins'  equations  produce  using  the  DELFIC  100  equa  1 -ac t  i  v  i t y 


groups  and  a  yield  of  1  megaton.  The  table  is  arranged  by 


decreasing  radius.  Again,  the  radius  is  the  mean  radius  of  the 


particles  in  the  group  whose  total  activity  comprise  one  percent 


of  the  total  activity  produced  by  the  1  megaton  explosion.  A 


complete  table  is  presented  in  Appendix  A. 


Table  II 

Hopkins'  Stabilized  Cloud  Data 


GROUP  RADIUS 

(pm) 


(km)  (km) 


GROUP  RADIUS 
(pm) 


(km)  (km) 


1  1917.0 

10  297.0 


148.0 

89.5 

57.5 
37.7 


0.0* 

7.608 

10.648 

11.841 

12.494 

12.898 


0.206 

1.690 

1.826 

1.880 

1.909 

1.927 


30.5 

24.5 
15.4 

8.89 

4.19 

.473 


13.044  1.934 

13.167  1.939 


13.352  1.948 

13.485  1.954 


13.581  1.958 

13.657  1.961 


•  iiii  ii 

Radii  for  groups  1-3  produce  unrealistic  (minus)  hc. 


It  is  clear  that  the  Hopkins  equations  produce  very 


similar  results  to  the  WSEG-10  equations.  This  is  to  be 


expected.  However,  by  inspection  of  the  55th  group,  above,  it 


appears  as  if  the  Hopkins  model  is  predicting  a  cloud  wherein 


the  activity  is  shifted  somewhat  lower  than  the  WSEG-10  model. 


Determining  the  Best  Cloud  Model 


To  determine  if  there  are  other  significant  differences, 


Xv'  s-'.y  v 


the  two  cloud  models  will  be  graphically  compared.  To  do  this, 
the  WSEG-10  model  will  be  displayed  as  a  vertical  gaussian  with 
the  median  at  an  hc  of  13.043  km  and  a  oz  of  2.348  km.  Figure  4 
displays  the  WSEG-10  model  of  the  activity  of  a  1  Mt  device 
versus  altitude. 

The  Hopkins  model  will  be  displayed  as  the  sum  of  all  100 
distributions  as  described  by  equation  (8).  Both  distributions 
will  be  expanded  in  the  horizontal  direction  to  assist  in  seeing 
the  correlation.  The  activity  will  be  spread  vertically  for  each 
group  using  Conners'  development  of  Hopkins'  equations.  Conners 
found  that, 

Azmi=ID+2rml3D 

ozio.25Azm  Cm)  (12) 

where 

Azm ^  is  the  predicted  thickness  of  the  ittl  group, 

ID  is  the  intercept  of  the  radius  in  meters, 

S|}  i3  the  slope  of  the  altitude  per  meter, 

Ojj1  i3  the  characteristic  width  of  the  ith  8rouP» 
and  both  Sq  and  Iq  are  given  by 

Sd  =  7-EXP(  1 . 78999-0. 0482491n(Y)«-0.0230248ln(Y)2 

-0.002259651n(Y)3+0.0001615191n(Y)4)  (13) 

Id  =  EXP (7.0351 8+0 .  1 589 1  41n ( Y ) +0 . 08375391n ( Y ) 2 

-0. 01554641n(Y)3+0. 000862 1 031n (Y) 4  )  (14) 

where  Y  is  the  yield  in  kilotons  (8:19). 


<M>t>  scaoxo/iiyi 


Figure  4.  Plot  of  Activity  vs  Altitude  Using  WSEG-10's  Model 
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By  treating  the  constituent  parts  separately,  Hopkins' 
model  predicts  a  vertical  activity  distribution  which  more 
closely  approximates  the  overall  shape  of  the  activity-plot 
derived  from  NCTS  rocket  data,  presented  by  WSEG-1G,  than  does 
the  single  Gaussian  approximation  using  the  WSEG-10  model. 
Therefore,  it  is  believed  that  the  Hopkins  model  for  vertical 
activity  is  better  than  the  WSEG-10  model.  WSEG-10's  model  for 
lateral  spread,  for  the  reasons  stated  previously,  remains  the 
best  model  for  the  other  spatial  distributions. 

It  is  cumbersome  to  deal  with  the  100  different 
characteristic  widths  associated  with  Conners'  technique.  To 
simplify  the  model,  the  plot  of  activity  versus  altitude  from 
Hopkins  (Figure  5)  is  fit  to  the  WSEG-10  plot  (Figure  4).  This 
is  done  byholding  the  area  under  the  Hopkins  plot  constant  and 
increasing  oz  until  the  three  inflection  points  on  both  curves 
coincide.  (See  Figure  6).  These  three  inflection  points  define 
°z(B)  for  the  bottom  half  of  the  plot,  the  peak  or  maximum 
horizontal  extent, and  the  oz(x)  *"or  tlie  toP  of  the  plot. 
The  values  of  the  two  oz  are  not  expected  to  be  the  same  since 
the  plot  is  clearly  skewed,  but  they  should  be  similar.  Also, 
because  the  simplification  will  be  defined  by  only  one  oz, 
fitting  the  two  plots  together  produces  an  area  of  underlap  in 
the  lower  altitudes  and  an  area  of  overlap  in  the  upper 
altitudes.  However,  the  central  portion  of  the  Hopkins  plot  can 
be  brought  into  coincidence  with  the  WSEG-10  plot.  The  overall 
oz  that  produces  this  coincidence  is  the  desired  oz  that  best 
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describes  the  fitted,  simplified  Hopkins  model  for  vertical 
activity.  The  oz  that  produces  the  degree  of  fit  evident  in 
Figure  6  is  2.28  km. 

This  completes  the  model  that  establishes  the  initial 
conditions  for  the  nuclear  cloud.  The  fallout  particles  fall 
to  the  ground  from  their  initial  spatial  positions  as  defined  by 
two  WSEG-10  lateral  Gaussian  distributions  and  one  Hopkins 
vertical  distribution  with  an  ammended  o2  of  2.28  km. 


1 


VI .  Phase  II  -  Fallout  Transport  in  the  Atmosphere 


A  single  object  falling  through  the  atmosphere  is  rela¬ 
tively  easy  to  follow  without  the  use  of  computers.  Wher^large 
numbers  of  objects  are  being  followed,  a  computer  is  obviously 
more  appropriate.  As  the  study  of  fallout  involves  many 
particles,  it  is  only  natural  that  several  well  defined  computer 
code  types  have  evolved.  The  two  most  common  types  are  the 
"Smear"  code  and  the  "Disk  Tosser."  Of  these  two  types,  the  disk 
tosser  is  more  adapted  to  following  the  specific  path  taken  by  a 
particle  as  it  falls. 

A  disk  tosser  begins  with  a  set  of  discrete  radii 
particles  representing  an  equal  mass  subdivision  of  the  particle 
size  distribution.  Each  particle  has  a  vertical  coordinate  which 
uniquely  defines  its  location.  The  cloud  model,  as  parameterized 
in  the  preceding  chapter,  is  this  set  of  coordinates.  Each  of 
the  100  particles,  then,  is  represented  by  an  initial  starting 
altitude.  A  disk  is  formed  by  taking  the  characteristic  width  of 
the  cloud  center  at  stabilization  and  using  it  as  the  radius  of 
the  disk.  The  amount  of  activity  contained  in  the  disk  depends 
on  its  distance  from  the  median  of  the  vertical  distribution 


j 


used  to  define  the  vertical  activity.  Because  the  effects  of 


l 


wind  are  ignored  prior  to  stabilization,  there  is  no  initial  J 

\ 

lateral  translation  of  the  disk  from  the  vertical  axis.  The  disk  J 

i 

tosser  then  applies  whatever  wind  conditions  are  considered  ! 


appropriate,  "tosses"  each  disk  into  the  wind  field,  and  the 
operator  observes  where  the  disk  goes.  A  simple  analogy  might  be 
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to  consider  each  disk  as  an  infinitely  thin  frisbee  thrown  by  a 


mechanical  arm  at  a  predetermined  height.  Clearly,  the  design  of 
the  wind  field  has  a  large  impact  on  the  path  of  the  disk.  This 
is  now  examined  in  more  depth. 

Modeling  the  Winds 

The  "X"  Wind.  By  definition,  the  direction  of  the  X  wind 
is  the  direction  of  the  wind  at  the  cloud  center  stabilization 
altitude.  The  direction  at  hc  is  deliberately  chosen  as  the 
downwind  direction  because  the  disk  that  starts  its  fall  from  hc 
will  have  the  maximum  amount  of  activity.  By  observing  the 
locations  of  all  the  disks  that  started  from  hc.  a  line  of 
maximum  activity  can  be  traced  on  the  ground.  This  line  is 
referred  to  as  the  "Hotline." 

In  the  simplest  case,  where  there  are  no  crosswind 
components,  the  hotline  and  X  axis  are  identical.  If  a  cross- 
wind  component  other  than  zero  is  chosen,  the  hotline  and  the  X 
axi3  will  diverge  from  a  common  origin,  i.e.  ground  zero. 

The  wind  speed  at  hc  is  only  critical  in  that  it  must  be 
non-zero.  Clearly,  a  zero  downwind  velocity  would  cause  no 
translation  in  the  downwind  direction.  Likewise,  any  positive 
non-zero  velocity  will  translate  the  fallout  in  the  downwind 
direction  according  to  the  following  equation, 


where 


ta  is  the  time  of  arrival, 

vx  is  the  downwind  wind  velocity,  (previously  W) 
x  is  the  downwind  translation. 

Again,  since  the  spatial  variables  are  independent  of  one 
another  and  the  downwind  translation  as  shown  by  equation  (15) 
is  not  a  function  of  Y,  the  downwind  wind  velocity  can  be 
ignored . 

The  "2"  Wind.  The  stratified  atmosphere  assumption 
allows  uniform  treatment  of  the  effects  of  air  density  on  fall 
velocity  with  altitude.  Inclusion  of  a  wind  with  a  Z  component 
alters  the  composition  of  the  atmosphere  as  related  by  U.S. 
Standard  Atmosphere.  WSEG-10  chooses  to  ignore  vertical  winds  by 
concluding  that  the  transport  is  "...  primarily  horizontal  not 
vertical."  (5:3)  More  specifically,  defining  a  non-zero  vertical 
wind  causes  the  otherwise  independent  spatial  variables  to  lose 
their  independence.  Panof3ky  and  Dutton  suggest  that  even  in 
the  more  complex  applications,  where  precise  expressions  of 
three  dimensional  wind  fields  are  desired,  the  three  variables 
are  treated  independently  (13:230-235).  Therefore,  vertical 
winds  are  not  treated  by  this  model. 

The  "Y"  Wind.  The  "Y"  wind,  or  crosswind,  acts  on  the 

particlesin  a  direction  normal  to  the  downwind  direction.  Since 

/ 

the  direction  of  the  wind  at  hc  iS  defined  as  the  downwind 
direction,  the  magnitude  of  the  normal,  or  Y  component  of  the 
wind  at  h_  mu3t  be  zero.  The  magnitude  of  the  crosswind  is  known 


to  change  with  altitude.  A  typical  value  of  the  change  lies 
between  1.0  and  10.0  km/hr-km  of  altitude.  It  is  this  change  in 
magnitude  with  altitude  that  is  defined  as  horizontal  shear.  Two 


examples  are  used,  next,  to  assist  in  refining  the  range  of 
shear  values  to  a  single  value  to  be  used  in  the  model  as  an 
average  shear. 


Actual  Observations.  On  5  July  1957,  in  the  Nevada 
desert,  a  74  Kt  device  was  exploded  as  the  final  portion  of 
Operation  Plumbbob  Hood.  Table  III  represents  a  portion  of  the 
wind  data  recorded  at  H-hour.  Figure  7  is  a  top-down  view  of  the 
data  from  Table  III  with  all  the  wind  vectors  acting  at,  or 
directly  above,  the  origin.  By  averaging  the  change  in  wind 
speed  from  the  surface  to  cloud  top,  one  can  see  that  the 
average  shear  is  6.32  km/km-hr.  Initially,  this  appears  to  be 
a  valid  shear  value  since  it  lies  within  the  range  of  typical 
shear  values  shown  above.  However,  since  the  cloud  is  the  result 
of  only  a  74  Kt  yield  device,  the  cloud  top  does  not  project 
beyond  the  tropopause.  The  shear  values  in  Table  III  naturally 
reflect  the  greater  turbulence  of  the  troposphere. 

The  second  example,  provided  by  the  Air  Force  Operational 
Test  and  Evaluation  Center  (AFOTEC),  helps  to  clarify  whether 
the  average  shear  developed  from  the  first  example  is  too  high 


or  not. 


AFOTEC  Wind  Projections.  AFOTEC  has  a  data  file  of 


observed  winds  at  various  pressure  heights  for  a  5  by  2.5  degree 
grid  that  spans  the  globe.  Accessing  this  data  file  for  the  wind 


Table  III 

Nevada  Wind  Data  for  5  July  1957 
Site  elevation  -  4230  ft  :  37°  08'  N,  1160  02'  W 


Altitude^!  Angle  Speed  Del  Altitude 

eg)  (mph)  (mph)  t  |(k-ft! 


II 


ground 

■59 

0 

6.0 

mi 

5 

8.0 

Bn 1 

9 

10.0 

■8 a 

12 

12.0 

150 

14 

14.0 

270 

23 

16.0 

180 

33 

18.0 

180 

24 

20.0 

180 

12 

22.0 

220 

9 

24.0 

230 

13 

Del  V 
(mph ) 


14  30.0 

15  32.0 

16  34.0 

17  36.0 

18  38.0 

19  40.0 

20  42.0 

21  44.0 

22  46.0 


notes:  1.  All  altitudes  are  with  respect  to  Mean  Sea  Level. 

2.  Cloud  top  -  48000  ft:  Cloud  bottom  -  35000  ft. 

3.  Tropopause  height  measured  at  53149  ft. 

4.  Change  in  speed  to  next  higher  altitude. 


(14:280-283) 


Scale:  1"  =  5  mph 

note:  Number  at  vector  head  refers  to  item  ( 0)  above 


Figure  7.  Observed  Nevada  Shear  above  Ground  Zero 


data  that  most  closely  matches  the  date  and  geographical 
location  of  the  Plumbbob  test  produces  the  data  shown  in  Table 
IV  for  16  July.  (The  data  represents  five  year  averages  of 
observations  and,  thus,  is  not  tied  to  a  particular  year.)  This 
data  is  depicted  in  Figure  8  in  the  same  manner  that  was  used 
for  Figure  7.  The  average  shear  is  found  to  be  0.18  km/hr-km. 
The  last  two  pressure  levels  shown  in  Table  IV  reflect  altitudes 
that  are  above  the  tropopause  as  measured  on  5  July  1957.  The 
last  pressure  level  (70  mb)  is  clearly  within  the  stratosphere. 
Here,  the  magnitude  of  the  wind  has  fallen  sharply  from  the 
levels  it  was  at  in  the  troposphere.  This  drop  in  magnitude  is  a 
contributor  to  the  lower  average  shear  value. 

Although  no  direct  comparison  of  these  two  examples  can 
be  made,  they  can  be  U3ed  to  help  bound  the  appropriate  value  of 
average  shear  for  the  model.  The  observed  data  covers  the  winds 
in  the  troposphere  in  detail  for  one  specific  day.  The  predicted 
data  from  AFOTEC  is  an  average  of  five  days  and  extends  to  a 
greater  height  than  the  observed  data.  The  "best"  value  lies 
between  the  two. 

The  WSEG-10  authors  use  a  shear  that  ranges  from  0.1  to 
0.4  knot3/1000  ft.  A  comment  is  made  that  these  shear  values  are 
more  common  at  50  k-ft  and,  that  at  a  lower  altitude  of  30  k-ft, 
higher  shear  values  are  found  (5:47).  Both  for  ease  of 
comparison  with  WSEG-10  results  and  as  a  median  value  of  the 
shear  values  discussed  above,  this  study  uses  1.0  km/km-hr  as 


the  effective  shear  for  all  altitudes.  A  200  meter  increment  of 


Table  IV 


altitude  and  0.2  km/200  m-hr  allows  a  smoother  transition  of 


wind  velocity  with  altitude  than  does  the  1.0  km/km-hr  shear 
value. 

Finally,  both  the  actual  shear  data  and  the  predicted 
shear  data  demonstrate  the  large  variability  in  direction  at 
different  altitudes.  A  particle  subject  to  this  variation  would 
find  itself  moving  back  and  forth  as  it  fell  to  the  ground.  The 
crosswind  distancet r a v e 1 1  ed  is,  therefore,  less  than  if  the 
crosswind  were  from  only  one  direction.  The  maximum  total 
influence  of  the  crosswind  is  achieved  by  allowing  it  to  act  on 
the  fallout,  in  the  same  direction,  at  all  altitudes.  To  model 
this  maximum  extent,  and  yet  still  constrain  the  value  of  shear 
at  hc  to  0.0,  a  convention  is  established  that  defines  the  shear 
above  hc  as  acting  in  the  negative  direction  and  below  hc  as 
acting  in  the  positive  direction  with  respect  to  the  Y  axis. 
Figure  9  illustrates  how  this  convention  operates  in  two 
dimensions  to  produce  a  discrete  yet  continuous  crosswind 
magnitude  in  a  case  where  h^  is  13.0  km. 

The  la3t  group  of  parameters  with  which  the  model  must 
contend  describe  the  downward  motion  of  the  particles.  These 
forces  must  be  balanced  to  produce  a  realistic  time  and  speed  of 


Altitude 


Figure  9.  Maximum  Shear  Effect  Model 

The  Fall  Mechanics 

Each  particle  of  fallout  must  remain  in  the  layers  of  the 
stratified  atmosphere  for  a  finite  amount  of  time  in  order  for 
the  crosswind  to  have  an  effect.  Considering  a  particle  at  a 
specific  time  within  its  total  fall  time,  leads  to  a  static 
condition  where  all  the  forces  acting  on  the  particle  are  in 
balance.  For  the  moment,  the  lateral  forces  exerted  by  the 
crosswind  will  be  ignored.  This  leaves  only  the  vertical  forces. 


It  is  possible  to  approximate  the  forces  by  beginning 


with  a  term  for  the  buoyant  force  as 


«  =  PDCd» 


where 


F4  is  the  upward  force  (N), 

Pp  is  the  dynamic  pressure  (kg/m-sec2), 

Cd  is  the  drag  coefficient  (dimensionless), 
o  is  the  cross  sectional  area  (m2). 


And,  the  downward  force  is  simply 


F*=Wg 


where 


F+  is  the  downward  force  (N), 

W  is  the  mass  of  the  particle  (kg), 
g  is  the  gravitational  constant  (m/sec2). 


The  assumption  that  all  of  the  fallout  particles  are 
spheres  allows  some  simplification  of  the  above  equations  as 
follows:  the  dynamic  pressure  on  a  sphere  is  Pp  =  (1/2)pavz2; 
the  cross  sectional  area  of  a  sphere  is  o  =  it  r g;  and  the  mass 
is  W  =  (4/3)irr3p For  the  forces  to  be  balanced,  F4  must  equal 
Fi .  Setting  them  equal  and  making  the  substitution,  one  gets 


( 1  /2)pav  2Cd* r2=(  4/3  )it  r3pfg 


where 


pa  is  the  density  of  the  air  (kg/m3), 
vz  is  the  velocity  of  the  particle  (m/sec), 
r  is  the  radius  of  the  particle  (m), 

Pf  is  the  density  of  the  fallout  (2600  kg/m3). 

From  equation  (18)  one  could  solve  for  the  velocity  by 
rearranging  the  terms  and  substituting  the  correct  values  for  pa 
and  Cj.  Unfortunately,  is  a  derived  value  based  on  the  fall 
velocity  of  the  particle  and  the  particle's  shape.  Fortunately, 
by  combining  the  work  of  Davies  and  McDonald,  as  presented  by 
Bridgman  and  Bigelow  (16:212),  the  drag  coefficient  can  be 
treated  in  such  a  manner  as  to  eliminate  the  need  to  solve  for 
Cd  explicitly.  Davies  shows  two  relationships  (equations  (19) 
and  (20))  between  Cd  and  Re  (Reynolds  number)  and,  Davies  shows 
a  relationship  (equation  (21))  between  Re  and  vz  (17:463).  These 
relationships  are, 

Re=Re2cd/24-2.3363El0“4(Re2Cd)2 

+2.0154El0-6(Re2Cd)3-6.9105El0-9(Re2cd)4  (19) 

for  Re2cd< 120 

log(Re)  =  -1 .29536t-0.986log  (Re2Cd) -0.046677  [log  (Re2Cd)32 

+  0.001  1235[log(Re2Cd)  ]3  for  Re2Cd>  =  120  (20) 

where  Reynolds  number  is  defined  as; 

Re= ( 2vzr )/ v  (21) 
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where  v  is  the  kinematic  viscosity  of  the  air  (m2/sec).  However, 
since  v=  n /pa  the  Reynolds  number  can  be  restated  as 

Re=(2vzpar)/n  (22) 

where  n  is  the  dynamic  viscosity  of  the  air  (kg/sec-m). 

Solving  equation  (22)  for  vz  yields 

vz=Ren/(2par )  (23) 

and  substituting  this  into  equation  (19)  McDonald  obtains  (17: 
463) 


Re^Cd= ( 32paPf gr ^ ) / ( 3n^ )  (24) 

Now,  the  quantity  Re^c^  is  defined  by  parameters  which 
are  known.  Re^Cj  can  be  used  in  equations  (19)  and  (20)  to 
determine  Re.  With  Re  defined,  equation  (23)  becomes  an 
expression  for  the  instantaneous  velocity  at  a  given  altitude. 
A  correction  factor  for  high  altitude  "slip"  is  suggested  by 
Davies.  Here  the  correction  factor  used  by  DELFIC,  as  reported 
by  Bridgman  and  Bigelow  (16:212),  is  used.  It  is 

CF= 1+ 1 . 165ElO”7/(rpa)  (25) 

The  final  expression  for  the  balanced  forces  on  the 
fallout  as  it  falls  is,  then. 


vz  =  (CF)RenA2pfr ) 


(m/sec ) 


(26) 


An  Expression  for  Fall  Time 

Equation  (26)  expresses  the  velocity  at  a  specific 
altitude.  In  order  to  arrive  at  an  equation  that  expresses  the 
time  of  fall,  a  few  additional  steps  must  be  taken.  Since,  in 
general , 


(velocity)(time)=distance 

if  the  distance  is  the  total  height  from  which  the  particle  must 
fall  and  velocity  is  the  velocity  in  the  z  direction,  then 

vzt=z  (27) 

and 

t=Az/vz  (28) 

and,  if  rather  than  take  the  entire  height  of  fall,  one  takes  z 
(in  equation  (28))  to  be  an  increment  of  altitude,  and  the 
velocity  is  assumed  constant  over  this  increment  of  altitude, 
then , 


where 


t=  A z/vz  (sec)  (29) 


t  is  the  amount  of  time  to  fall  through  Az, 
Az  =  200  m, 

vz  is  defined  by  equation  (26)  (m/sec). 


If  equation  (29)  provides  the  time  of  fall  through  one  increment 


of  altitude,  then  adding  the  times  from  all  the  increments  gives 
the  total  time  of  fall.  The  total  time  of  fall  is,  clearly,  the 
time  of  arrival  of  the  particle  on  the  ground,  and  it  is 
referred  to  as  ta.  stated  another  way, 

n 

ta=  ^](Az/vzi)  (sec)  (30) 
i  =  1 

By  combining  equation  (29)  and  the  shear  force  at  each 
altitude,  the  lateral  translation  of  the  particle  is  determined. 
Ultimately,  as  the  particle  is  transported  through  the  lowest 
increment  of  altitude  and  is  translated  by  the  final  shear 
force,  the  particle  arrives  at  its  "grounded  location." 

Figure  10  is  a  graphical  representation  of  these  balanced 
forces  with  a  time  history  of  a  hypothetical  particle.  So  that 
the  mechanism  is  more  readily  visible,  the  particle  is  only 
lofted  to  a  height  of  6  km.  Cloud  center  stabilization  altitude 
is  chosen  to  be  4  km  so  that  the  effect  of  the  shear  convention 
above  hc,  in  the  negative  direction,  is  more  apparent.  The 
horizontal  shear  is  1  km/km-hr.  For  the  purpose  of  this 
illustration  only,  the  atmosphere  is  considered  to  be 
homogeneous  so  that  the  fall  velocity,  vz,  becomes  constant.  The 
showing  of  the  "fall  distance"  arrow  as  a  force  vector  is  merely 
an  artifice  to  assist  visualization.  The  "effective  wind"  is  any 
arbitrary  non-zero  wind  with  units  of  km/hr.  A  "track"  in  the 
X-Y  plane  is  the  ground  projection  of  the  "path  travelled."  This 
is  representative  of  what  the  computer  code  is  actually  storing 
as  the  lateral  distance  of  the  particle  from  the  downwind  axis. 


Fiaure  10.  The  Tran. - 
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With  a  set  of  values  placing  a  particle  in  the  cloud 
model  as  an  initial  condition,  the  transport  model  operates  to 
determine  that  particle's  grounded  location.  The  third  phase  of 
pattern  formation  requires  the  grounded  location  of  all  the 
particles  to  be  determined.  However,  before  this  phase  is 
discussed,  some  intermediate  results  from  the  examination  of 
just  one  or  two  particles  will  be  discussed. 

Intermediate  Results 

The  interface  between  the  cloud  model  and  the  transport 
model  is  accomplished  by  use  of  Colarco's  equation  and  coeffi¬ 
cients  (18:11,67-71).  Colarco  shows  that  by  replacing  the  verti¬ 
cal  distribution  by  a  Dirac  delta  function,  the  following 
equation  predicts  the  particle  size  that  arrives  at  tg. 

6 

r  =£  ^ita^”^+c7ta”’ ^  (urn)  (31) 
i=  1 

where  and  Cy  are  the  coefficients  from  Colarco's  work  andta 
is  in  hours.  When  the  coefficients  for  an  altitude  of  13  km  are 
used,  the  resultant  radii  define  those  particles  that  originated 
at  hc.  This  is  particularly  useful  for  determining  the  location 
of  the  hotline.  Furthermore,  with  all  of  the  activity  of  the 
disk  confined  to  the  Dirac  function,  rather  than  distributed 
under  a  gaussian  curve,  the  deviation  of  the  grounded  location 
in  this  case  from  the  downwind  axis  is  the  translation  due  to 


In  the  process  of  applying  Colarco’s  equation  to 
determine  an  initial  set  of  data  points  defining  the  hotline  for 
a  stabilization  altitude  of  13*04  km  and  times  from  3  to  24 
hours,  it  was  decided  to  perform  a  secondary  calculation  using 
the  transport  model.  The  radius  predicted  by  equation  (31)  was 
entered  into  equation  (30).  Both  equation  (30)  and  (31)  should 
have  yielded  the  same  ta.  As  can  be  seen  in  Table  V,  the 
Colarco  predicted  radius  was  found  to  be  slightly  different 
than  the  radius  required  by  equation  (30)  to  achieve  the  same 
t_.  This  discrepancy  was  not  resolved.  Colarco  used  a 

a 

polynomial  fit  routine  to  determine  the  coefficients.  Because  of 
this  method,  it  is  believed  that  some  degree  of  error  was 
introduced  to  the  coefficients.  Therefore,  equation  (30)  was 
thought  to  be  more  accurate.  Colarco's  method  proved  to  be  a 
quick  way  of  approximating  the  correct  radius.  The  Colarco 
prediction  was  then  incrementally  adjusted  until  the 
radius,  after  substituting  into  equation  (30),  produced  the 
desired  ta* 

As  stated  previously,  the  downwind  wind  velocity  is  not 
critical.  Merely  as  a  point  of  reference,  Table  V  shows  the 
downwind  translation  using  a  vx  of  50  km/h r. 

Just  as  the  particle  originating  from  hc  defines  the 
hotline,  so  also  does  the  hotline  define  the  center  of  the 
lateral  distribution  on  the  ground  in  the  Y  direction  when  the 
activity  is  not  constrained  by  a  Dirac  function.  This  follows 
from  the  assumption  that  the  vertical  distribution,  which  is 
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Table  V 

1  Mt  Hotline  Locations 


II 


_  1  2  Y 

rc  r  m  A 
(um)  (pm)  (km) 


3 

75.00 

73.80 

150 

4.38 

4 

61.43 

59.78 

200 

5.83 

6 

47.  12 

45.89 

300 

8.64 

8 

39.35 

38.54 

400 

11.52 

10 

34.33 

33.72 

500 

14.37 

12 

30.75 

30.32 

600 

17.20 

14 

28.04 

27.78 

700 

20.03 

16 

25.91 

25.81 

800 

22.85 

18 

24.16 

24.20 

900 

25.68 

20 

22.70 

21.73 

1000 

28.51 

22 

21.47 

21.22 

1100 

31.35 

24 

20.39 

20.76 

1200 

34.16 

13.15 

17.48 

25.93 

34.57 

43.10 

51.60 

60.09 

68.56 

77.03 

85.52 

94.05 


17.54 

21.92 

23.30 

29.13 

34.57 

43.21 

46.09 

57.61 

57.46 

71.83 

68.80 

85.99 

80.12 

100.15 

91.41 

114.26 

102.71 

128.38 

114.02 

142.53 

125.40 

156.75 

136.63 

170.79 

77 

65 

28 


23.04 

28.73 

34.40 

40.06 

45.71 

51.35 

57.01 

62.70 

68.32 


notes:  1.  rc  is  the  Collarco  predicted  radius, 

2.  rm  is  the  adjusted  (correct)  radius, 

3.  All  shear  have  units  of  kra/km-hr. 


approximated  by  the  gaussian  distribution,  is  translated  to  a 
gaussian  distribution  in  Y.  Subsequently,  a  particle  which  is 
defined  as  being  a  distance  oz  from  hc  will,  after  being 
transported  to  the  ground,  be  a  distance  defined  by  oy  from  the 
hotline.  The  importance  of  the  shape  and  location  of  the 
hotline,  then,  is  sufficient  to  warrant  a  closer  look. 

Figure  11  shows  the  predicted  locations  of  the  hotlines 
formed  by  shears  of  0.2,  0.6,  and  1.0  according  to  the  data 
listed  in  Table  V.  Clearly,  in  the  regime  from  3  to  24  hours, 
when  only  the  translation  due  to  shear  is  considered,  the 
hotlines  project  as  straight  lines. 

Following  the  same  procedure  that  was  used  to  determine 


the  hotlines,  the  grounded  locations  of  the  particles  that 


Table  VII 

Oy  Values  for  5  Shear 

fca 

.2 

.4 

.6 

.8 

1.0 

(hr ) 

(km) 

(km) 

(km) 

(km) 

(km) 

■  -L  1  — U  —  I   - □ 

3 

.580027 

1 .16005 

1.74007 

2.32001 

2.90013 

4 

.783387 

1.57678 

2.36516 

3.15354 

3.94193 

6 

1 .23034 

2.46067 

3.69101 

4.92135 

6.15168 

8 

1.65997 

3.31993 

4.97989 

6.63986 

8.29983 

10 

2.07494 

4.14989 

6.22484 

8.2998 

10.3747 

12 

2.51727 

5.03454 

7.55181 

10.0691 

12.5864 

14 

2.96347 

5.92694 

8.89041 

1  1.8539 

14.8173 

16 

3.41088 

6.82175 

10.2326 

13.6435 

17.0544 

18 

3.84931 

7.69863 

1  1.5479 

15.3972 

19.2466 

20 

4.29521 

8.59043 

12.8856 

17.1808 

21.4760 

22 

4.75827 

9.51655 

14.2748 

19.0331 

23.7913 

24 

5.18831 

10.3766 

15.565 

20.7533 

25.9416 

note:  All  3hear  in  units  of  km/km-hr 


Figure  12.  Plot  of  Oy  versus  Time 
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computed  difference  between  the  locations  for  the  five  values  of 
shear  used  for  Tables  V  and  VI.  Figure  12  is  a  plot  of  these  oy 

Before  attempting  to  look  at  the  patterns  formed  by  all 
the  particles,  the  vertical  increment  of  fall  distance  will  be 
analyzed  in  an  effort  to  ensure  that  it  is  the  most  appropriate 
one  to  use  for  the  completed  computer  model. 

Variations  of  the  A  z  Parameter 

The  Dirac  function  substitution  provides  a  quick  method 
of  testing  one  of  the  primary  assumptions.  The  use  of  200  meters 
as  the  "best"  increment  of  vertical  distance  is  used  both  for 
the  shear  increment  and  the  fall  distance  over  which  the  fall 
velocity  is  held  constant.  A  comparison  of  this  parameter 
establishes  the  validity  of  this  assumption. 

Table  VIII  shows  the  results  of  varying  the  vertical 
increment,  Az,  from  1000  m  to  50  m.  50  m  is  considered  to  be  the 
most  accurate  increment  available  because  the  parameters 
defining  vz.  from  U.S.  Standard  Atmosphere,  are  not  given  in  any 
finer  detail.  Therefore,  the  time  of  fall  that  is  computed  using 
Az  =  50  m,  is  assigned  a  degree  of  accuracy  of  100%.  The  other 
Az  are  compared  to  the  50  m  fall  time  and  assigned  a  relative 
degree  of  accuracy  accordingly. 

Deciding  which  Az  to  use  is  largely  subjective.  However, 
Table  VIII  demonstrates  that  the  primary  assumption,  that  200  m 
is  sufficiently  accurate  to  predict  valid  results,  is  acorrect 
one.  The  most  accurate  result  is  achievable  with  Az  =  50  m,  but 


VII.  Phase  III  -  The  Grounded  Fallout  Pattern 


The  grounded  fallout  pattern  is  a  result  of  plotting  the 
locations  of  all  the  disks  that  comprise  the  nuclear  cloud.  The 
finer  the  mesh  of  disks  in  the  cloud  is  made,  the  more  accurate 
the  pattern  on  the  ground  will  be.  The  DELFIC  program,  for 
instance,  uses  100  disks  to  describe  the  vertical  distribution 
of  the  activity.  Each  of  these  disks  is  then  further  subdivided 
into  100  additional  disks  with  a  representative  particle  size 
describing  0.01S  of  the  original  activity  of  the  initial  disk. 
This  means  that  the  DELFIC  code  actually  follows  10,000  separate 
disks.  This  mesh  is  considered  to  be  fine  enough  to  predict 
relatively  accurate  fallout  patterns.  However,  even  mainframe 
computers  need  considerable  time  to  handle  the  data  required  for 
such  a  fine  mesh.  But,  when  only  the  relative  amplitude  across 
the  pattern  at  any  downwind  location  is  required,  as  it  is  for 
the  study  of  shear  growth  with  time,  only  a  sub  set  of  all 
possible  disks  must  be  analyzed.  This  simplifies  and  speeds  up 
the  process.  This  simplification  allows  the  program  to  run 
handily  on  even  a  microcomputer. 


Building  the  Pattern 


A  disk  tosser  code,  like  DELFIC,  does  not  permit  snap¬ 
shot  analysis  of  fallout  along  a  crosswind  axis.  When  the 


downwind  axis  is  considered  as  a  time  line,  as  in  Figure  10, then 
the  crosswind  axis  becomes  a  time-slice  of  the  developing 
fallout  pattern.  With  the  particles  constrained  to  fixed 
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starting  altitudes,  the  different  disks,  which  are  further 
defined  by  their  100  particle  radii,  fall  at  different 
velocities  as  determined  by  equation  (30).  This  means  that  no 
two  disks  arrive  on  the  ground  at  the  same  time.  Therefore,  the 
centers  of  the  grounded  disks  are  offset  in  both  X  and  Y.  This 
offset  makes  it  necessary  to  analyze  all  of  the  disks  which  have 
fallen  prior  to  the  chosen  time  of  analysis  so  that,  by  adding 
the  two  dimensional  circular  normal  functions  which  describe  the 
activity  in  the  grounded  disks,  the  correct  amplitude  at  any 
ground  location  can  be  found. 


A  two  step  sorting  process  is  used  to  overcome  the  need 
to  analyze  all  of  the  disks  mentioned  above.  Fir3t,  Colarco's 
equation  (31)  is  used  to  determine  the  approximate  particle  size 
which  arrives  at  the  arrival  time  of  interest.  Second,  the  radii 
predicted  by  equation  (31)  are  used  as  delimiters  in  an 
iterative  routine  which  processes  the  100  equal  activity  radii 
through  equation  (30)  with  altitude  parameters  established  by 
the  initial  stabilized  cloud  top  and  cloud  bottom.  This 
effectively  establishes  crude  starting  altitudes  for  the  full 
range  of  particles  that  arrive  close  to  the  desired  time  of 
arrival.  Table  IX  is  an  example  of  this  sorting  process  where 
cloud  top  was  set  at  18,000  m  and  cloud  bottom  was  set  at  1000 
m.  The  vertical  increment  is  1000  m. 

The  results  of  the  initial  sort  are,  understandably,  only 
crude  approximations.  The  process  of  sorting  is  continued  to 
whatever  degree  of  accuracy  is  desired.  As  the  sort  is 
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Table  IX 

Initial  Sort  of  Radii  to  Determine  Starting  Altitude 

Az  =  1000  m  :  Cloud  top  =  18000  m  :  Cloud  bottom  =  1000  m 
Input  arrival  time  for  equation  (31)  =  3  hrs 

Radius 

Altitude 

Radius 

ta 

Altitude 

(pm) 

(hrs ) 

(m) 

(urn) 

(hrs ) 

(m) 

1917 

.1956835 

18000 

71.4 

3.128213 

12000 

1064 

.2801108 

18000 

68.4 

3.111329 

11000 

782 

.3478209 

18000 

65.5 

3.076789 

10000 

629 

.4105415 

18000 

62.7 

3.020838 

9000 

529 

.4720622 

18000 

60.1 

3.2011 

9000 

457 

.5341228 

18000 

57.5 

3.114077 

8000 

403 

.5964031 

18000 

55.1 

3.304502 

8000 

361 

.6589445 

18000 

52.9 

3.163483 

7000 

326 

.7246893 

18000 

50.7 

3.368404 

7000 

297 

.792291  1 

18000 

48.6 

3.190648 

6000 

272 

.863637 

18000 

46.6 

3.419007 

6000 

251 

.9360802 

18000 

44.7 

3.183808 

5000 

232 

1.014603 

18000 

42.8 

3.421722 

5000 

216 

1.09313 

18000 

41.1 

3.090501 

4000 

202 

1.173736 

18000 

39.4 

3.311883 

4000 

189 

1.261139 

18000 

37.7 

3.558912 

4000 

177 

1.355295 

18000 

36.2 

3.083957 

3000 

167 

1.446057 

18000 

34.7 

3.305018 

3000 

157 

1.550671 

18000 

33.2 

3.55466 

3000 

148 

1.659398 

18000 

31.8 

3.818962 

3000 

140 

1.770279 

18000 

30.5 

3.  1  12994 

2000 

133 

1.880645 

18000 

29.2 

3.352316 

2000 

126 

2.006053 

18000 

28 

3.603631 

2000 

119 

2.149685 

18000 

26.8 

3.889906 

2000 

113 

2.290389 

18000 

25.6 

4.218101 

2000 

108 

2.422538 

18000 

24.5 

3.079243 

1000 

103 

2.570773 

18000 

23.5 

3.319904 

1000 

98.4 

2.723914 

18000 

22.4 

3.623523 

1000 

93.8 

2.895769 

18000 

21.4 

3.941959 

1000 

89.5 

3.076839 

18000 

20.5 

4.270033 

1000 

85.5 

3.017588 

16000 

19.5 

4.690053 

1000 

81.7 

3.06394 

15000 

18.7 

5.076425 

1000 

78.1 

3.09858 

14000 

17.8 

5.575694 

1000 

74.7 

3.  1  19044 

13000 

17 

3.081959 

4 

note:  The  17  pm  particle 

is  included  to  show  that, 

in  fact 

smaller  particles  require  less 

than  I 

the  cut* 

-off  altitude 

in  order  to  have  an  1 

arrival 

time  near 

3  hours. 
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continued,  those  radii  which  can  not  have  a  starting  altitude 
within  the  cloud  parameters  and  still  achieve  the  arrival  time 
desired,  are  rejected.  The  initial  sort  is  quickly  reduced  to  a 
group  of  radii  which  have  the  same  arrival  time.  Once  the  final 
sort  is  accomplished,  the  starting  altitude  is  used  with  the 
initial  distribution  data  from  equations  (11)  and  (12)  along 
with  the  standard  formula  for  a  gaussian  distribution  to 
determine  the  exact  activity  of  the  disk  represented  by  the 
particular  radius  at  that  altitude.  Since  the  total  activity  of 
any  one  of  the  100  radii  is  0.01%  of  the  total  cloud  activity, 
the  activity  present  in  one  of  the  disks  determined  by  this 
method  must  be  the  amplitude  of  the  gaussian  at  that  altitude. 
Table  X  shows  the  results  of  the  final  sort  using  an  altitude 
increment  of  1.0  m.  This  is  representative  of  the  data  that  is 
used  for  the  remainder  of  this  study. 


Table  X 

Final  Sort  of  Radii  for  ta=3hr :Determination  of  Activity-Ag1 


Group  Radius  Altitude  Activity | [Group  Radius  Altitude  Activity 
(um)  (m)  (Ci)  II  (pm)  (m)  (Ci) 


89.5 

85.5 

81.7 
78.1 

74.7 

71.4 

68.4 


8323 

1 .96E-6 

6863 

2.41E-5 

5553 

1 .28E-4 

4373 

3.60E-4 

3313 

6.24E-4 

2323 

7 . 55E-4 

1464 

6.92E-4 

- - J - 

65.5 

62.7 
60.  1 

57.5 
55.1 
52.9 

50.7 

10664 

9914 

9244 

8614 

8044 

7534 

7023 

5. 17E-4 
3.27E-4 
1.87E-4 
9 . 66E-5 
4.78E-5 
2.33E-5 
1.05E-5 

due  to  any  X  shear,  and  therefore  the  radius  of  the  grounded 
disks  remains  the  same  as  at  stabilization  height,  the  use  of  a 
non-zero,  positive  X-direction  effective  wind  guarantees  that, 
with  a  time  increment  of  1  hour  between  points  of  analysis, 
there  can  be  no  overlap  of  the  disks  in  the  X  direction.  This 
ensures  the  ability  to  analyze  the  crosswind  spread,  at  any  time 
t,  without  concern  that  portions  of  the  disks  that  landed  at  t-1 
are  overlapping  the  disks  at  time  t. 

Additionally,  as  the  time  of  arrival  increases  to  24 
hours,  the  quantity  of  disks  arriving  is  reduced.  This  is  a 
natural  result  of  the  initial  cloud  parameters  which  showed  that 
only  the  very  small  particles  would  stabilize  at  the  top  of  the 
cloud.  The  larger  particles  have  fallen  out  before  24  hours 
passes.  This  leaves  fewer  and  fewer  particles  airborn  as  time 
progresses.  An  example  of  this  reduction  is  demonstrated,  as 
shown  in  Table  XI,  in  .the  data  for  disks  arriving  at  24  hours. 


Table  XI 

Final  Sort  of  Radii  for  ta=24hr :Determination  of  Activity-Ag* 


Group  Radius  Altitude  Activity | [Group  Radius  Altitude  Activity 
(um)  (m)  (Ci)  ||  (pm)  (m)  (Ci) 


23.5 

22.4 

21.4 

20.5 


4.00E-4 
1 .04E-4 
1 .58E-3 
1.68E-3 


1 .33E-3 
8.99E-4 
4.91E-4 
2.54E-4 


When  the  particles  in  the  sorted  groups  are  operated  on 
by  the  transport  model,  they  form  a  line  ofpoints  in  the 


crosswind  direction  parallel  to  the  Y  axis.  The  single  particle 


which  defines  the  center  of  each  of  the  disks  lands  at  the 


determined  point.  The  point  is  expanded  in  the  X  and  Y 
directions  by  application  of  the  two  gaussian  distributions 
which  originally  defined  the  disk  at  stabilization.  (See  Figure 


The  lack  of  an  X  shear  leaves  the  gaussian  distribution 
defining  the  downwind  spread  of  activity  in  each  disk  unchanged. 
The  activity  in  the  disk  would  normally  have  to  be  defined  by  a 
circular  normal  probability  distribution,  where  the  amplitude  at 
any  point  is  (19:936), 


-1/2{ [(x-mx)/ox]2+[ (y_m  )/0  ]2} 

P ( x  , y  )  =  [  1  / ( 2 ifOy ox )  ]  e  (32) 

but, with  <7y  =ax  at  stabilization,  and  the  median  of  the  downwind 
distribution  constrained  to  be  on  the  crosswind  axis  by  the 
sorting  process,  where  m  =x,  equation  (32)  reduces  to 


„  -1/2[(y-mv)/o0r 

2n  .  1 


P(y)*[1/(2iro0*)]  e 


where 


y  is  any  point  along  the  radius  of  the  disk  (km), 
my  is  the  point  calculated  by  the  transport  model  (km). 


Also,  with  the  distribution  in  X  unchanged,  the  peak  of 
the  distribution  must  have  a  relative  amplitude  of  0.3989  (20: 
584).  This  allows  further  reduction  of  equation  (33)  to  a  single 
distribution  in  the  Y  direction.  The  amplitude  of  any  point 
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along  the  distribution  can  now  be  determined  by  integrating  all 
the  disks  along  the  crosswind  axis,  or, 

f+OO  0.39AQi  -1/2[(y-mv)/o0]2 

Act(y)=  /  —  — —  e  dy  (Ci)  (34) 

^  IX rJ  2lT'°0 

w h e r eAg i  is  the  initial  activity  of  the  disk  determined  by  the 
sorting  process. 

When  equation  (3*0  is  analyzed  from  minus  infinity  to 
plus  infinity,  it  produces  mathematically  complete  solutions. 
However,  realistically,  the  pattern  does  not  extend  to  infinity. 
The  model  uses  a  limit  of  3oo»  instead  of  infinity.  Since  this 
study  does  not  consider  an  actual  pattern,  the  3cjq  limit  is 
used  simply  to  speed  computation  and  assist  in  graphical 
presentation.  The  implementation  of  these  limits  is  demon¬ 
strated  on  a  specific  example. 

During  the  numerical  integration  along  the  crosswind 
axis,  the  model  invokes  some  fundamental  identities  of  the 
gaussian  distribution.  The  gaussian  distribution  is  centered 
about  the  median,  u  ,  with  681  of  the  total  area  under  the  curve 
between  y  -  o  and  y  +  o.  This  leaves  16X  of  the  total  area  in  the 
two  "tails"  (1:929).  Since  each  of  the  disks'  activity  is 
distributed  in  a  gaussian  manner,  su p e r p o s i t i o n  i  n g  the 
d  i  s  t  r  i  b  u  t  i  o  n  s  of  all  the  disks  also  forms  a  gaussian 
distribution.  Every  point  of  activity  determined  by  equation 
(31*)  is  the  amplitude  of  the  s  u  p  e  r  p  o  s  i  t  i  o  n  e  d  gaussian 
distribution.  As  the  activity  is  determined  along  the  axis,  the 


amplitude  is  used  in  a  trapezoidal  integration  routine  to 
determine  the  area  under  the  superposition.  Those  points  along 
the  crosswind  axis  which,  after  being  analyzed  to  determine  the 
amplitude  of  the  activity  there,  sweep  out  an  area  of  the 
superposition  equal  to  16%,  50%,  and  84%  are,  respectively,  the 

upwind  o y u ,  the  median  (or  hotline),  and  the  downwind  o yd* (Here, 
upwind  and  downwind  refer  to  the  crosswind  direction.)  Just  as 
there  was  a  difference  in  ozx  and  azg  in  the  model  of  the 
stabilized  cloud  (refer  to  Figure  6),  a  difference  in  oyu  and 
oyd  is  expected. 

As  with  the  two  previous  portions  of  the  complete  model, 
some  examples  of  the  crosswind  model  in  operation  are  examined 
next.  Figures  14,  15,  and  16  show  the  superpositioned  gausslan 
distribution  for  arrival  times  of  3,  12,  and  24  hours.  The  3hear 
value  is  0.2  km/hr/200  m.  The  left  limit  of  activity  is  obtained 
by  subtracting  3oq  from  the  grounded  location  of  the  upwind- 
most  particle  in  the  3  hour  arrival  group.  The  right  limit  of 
activity  is  obtained  by  adding  the  same  to  the  downwind-most 
particle  in  the  24  hour  arrival  group.  (  From  equation  (5),  a  1 
Mt  yield  results  in  o0  =  2.9376.)  The  amplitudes  are  uniformly 
increased  to  permit  viewing.  The  perspective  is  from  a  point 
along  the  Y  axi3.  Ground  zero  is  indicated  by  "GZ."  Finally, 
because  only  the  relative  amplitude  is  of  concern  in  this 
process,  no  amplitude  scales  are  included. 

The  data  used  for  these  plots  is  presented  in  Appendix  B. 

The  final  output  of  the  three  phased  model  is  composed  of 
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y;.' 

V.y  all  the  0y  calculated  by  the  crosswind  patterns.  It  is, 

therefore,  particularly  significant  to  understand  exactly  what 
values  the  crosswind  model  determines.  A  clear  interpretation 
leads  to  enhanced  understanding  of  the  final  output. 


Interpreting  the  Pattern 


Four  significant  items  are  demonstrated  by  the  crosswind 


activity  plots  that  affect  subsequent  analysis. 


Unlike  the  traditional  fallout  patterns  that  have  closed 


shapes,  the  analysis  based  on  the  mass  of  grounded  fallout  does 


not  appear  to  have  a  closed  configuration.  This  is  a  logical 


consequence  of  the  cloud  and  transport  models.  The  cloud  model 


would,  if  not  arbitrarily  constrained  to  some  cloud  top  limit. 


continue  to  predict  smaller  and  smaller  particles  at  higher  and 


higher  altitudes.  As  is  evident  by  Tables  V  and  VI,  the  smaller 


particles  will  be  carried  further  than  the  larger  particles.  In 


practice,  when  the  dose  rate  is  calculated,  the  particles 


landing  at  late  times  carry  very  little  of  the  initial  activity. 


The  fallout  patterns  that  predict  dose  rate  contours  thus  show 


negligible  dose  rates  for  late  times.  This  model,  on  the  other 


hand,  continues  to  carry  these  particles  to  ever  greater 


distances.  This  phenomenon  is  clearly  visible  in  the  three 


activity  plots  shown  in  Figures  14,  15,  and  16.  The  activity 


pattern  gets  increasingly  wider  with  time.  Because  the  total 


activity  is  conserved,  as  the  width  increases  the  amplitude  must 


decrease.  This  preserves  the  area  of  each  of  the  activity  plots 


mm 


at  unity.  However,  since  the  quantity  of  disks  landing  at  later 
times  is  reduced  from  that  landing  at  earlier  times,  the 
activity  plot  for  24  hours  can  be  seen  to  be  less  smooth  towards 
the  right  limit  of  the  pattern.  This  is  an  inherent  shortcoming 
of  the  discrete  mesh  used  by  all  disk  tossing  codes.  It  causes 
fluctuations  in  the  determination  of  activity  by  the  integration 
routine  for  times  greater  than  20  hours  which  must  be  considered 
when  analyzing  the  final  output. 

The  difference  in  o zj  and  o zg  is  expected  to  manifest 
itself  in  the  grounded  activity  plots.  An  unexpected  result  of 
analyzing  the  two  grounded  <jy  is  the  reversal  of  the  difference 
in  magnitude.  At  early  times  o yu  is  greater  thanoyq.  This 
difference  decreases  to  the  7  hour  mark  where  they  are  nearly 
equal. oyd  then  grows  larger  than  oyu.  (See  Table  XII,  below.) 
This  is  a  result  of  the  particle  activity  distribution.  The 
heavier  particles  at  the  bottom  of  the  cloud  fall  out  in  the 
early  activity  plots,  thus  creating  more  activity  to  the  left, 
or  upwind,  side  of  the  median  of  the  distribution.  At  later 
times  (ta  >  7  hours)  the  increasingly  smaller  particles  are 
carried  proportionately  further  from  the  hotline  on  the  downwind 
side.  This  stretches  the  pattern  out  in  the  downwind  direction 
and  increases  the  downwind  Oyj  correspondingly.  A  different 
particle-activity  distribution  would  exhibit  a  different  spread 

Of  Oy. 

Despite  the  aberrations  in  the  activity  plots  for  times 
greater  than  20  hours,  the  shape  of  the  activity  plots  verifies 


the  initial  assumption  that  the  crosswind  activity  is 
distributed  in  a  gaussian  manner.  Although  this  is  self  evident, 
had  a  different  curve  emerged  after  superpositioning,  subsequent 
analysis  would  have  had  to  have  been  modified  along  with  the 
initial  assumption. 


Table  XII 


Differences  in  oyu  and  ayd 
ilculated  for  a  shear  of  0.2  fcm/hr-km) 


(hrs ) 

Left  Particle 
Location  (km) 

Right  Particle 
Location  (km) 

i-  i  —r 

(&) 

Difference 

(km) 

3.4567 

4.6492 

5.3922 

6.4844 

7.2469 

8.3434 

9.7367 

11.4742 

12.0857 

12.8071 

13.7391 

14.7787 

15.9114 

17.3221 

18.7959 

18.035 

19.8703 

22.0551 

21.2974 

23.6951 

22.8833 

25.369 


5.5321 

7.1939 

9.2362 

10.8744 

13.0132 

14.4568 

16.4707 

18.0313 

20.4215 

22.1979 

23.9945 

24.8878 

26.9937 

29.0556 

31.1628 

32.0484 

34.5593 

36.8599 

37.8385 

40.5052 

41.491 

44.2221 


2.9785 

2.9973 

3.0375 

3.1410 

3.258 

3.33 

3.429 

3.529 

3.6495 

3.7710 

3.9059 

4.0005 

4.1625 

4.3064 

4.4190 

4,6440 

4.8150 

4.9680 

5.2694 

5.3820 

5.6024 

5.6565 


2.9065 

2.9263 

2.9805 

3.113 
3.2594 
3.344 
3.454 
3.560 
3.7085 
3.844 
3.9719 
4.0765 
4.2445 
4.4024 
4.526 
4.763 
4.953 

5.113 
5.4284 
5.554 
5.7944 
5.8715 


Finally,  the  transport  model  deposits  some  particles  for 
the  early  times  very  close  to  the  effective  downwind  axis.  This 
follows  from  the  assumption  that  torroidal  growth  can  be  ignored 
for  this  analysis.  Without  torroidal  growth  in  the  first  three 


VIII.  Comparisons  with  WSEG.  Norment  and  Bridgman 


WSEG-10  Growth  Predictions 


shear  using  the  definition  of  crosswind  shear  from  equation  (2). 
The  three  shear  values  are  0.2,  0.6,  and  1.0  km/hr-km. 


Figure  17.  WSEG-10  Lateral  Growth  Predictions 


Norment  Growth  Predictions 


Similarly,  Norment’s  predictions  can  be  plotted  for  the 
3ame  three  3hear  values  using  equation  (4).  Zj  and  zg  are 
computed  by  two  empirical-  formulas  shown  by  Norment  (6:19)  as, 


vvv::--: 


Figure  18.  Norment  Lateral  Growth  Predictions 


Bridgman.  C.  J.  Bridgman,  during  the  process  of  advising 
the  conduct  of  this  research  effort,  derived  a  third  equation 
for  the  growth  of  o3.  His  unpublished  derivation  is  introduced 
here  in  order  to  show  the  origin  of  the  final  equation  that  is 
used  as  a  "benchmark"  for  the  results  of  the  model  explained  in 
this  study. 

Bridgman's  derivation  begins  with  the  knowledge  that  the 
difference  between  the  grounded  location  of  a  particle 
originating  at  hc  and  the  grounded  location  of  a  particle 
originating  at  hc+oz  is  oy.  Then,  looking  at  the  particle  which 


2 


started  at  h 


c  * 


(37) 


Substituting  vz  =  SyLzg-ztt)]  into  equation  (37)  and  making  a 
change  of  integration  limits  using  dt  =  -(dz/vz)  yields. 


Y1 


(zq-z) (-dz/vz) 


(38) 


He  then  makes  the  assumption  that  the  change  in  the  fall 
velocity  over  the  range  of  the  particle’s  fall  is  small  compared 
to  the  initial  velocity.  This  allows  replacement  of  the  vz  in 
equation  (38)  with  an  average  fall  velocity,  <vz>.  Now, 
integrating,  and  combining  terms  yields, 

Y i  =  (Syz02)/(2<vz>)  (39) 

Letting  <vz>  =  zg/t  in  equation  (39)  then  produces 

Y-i  =  Sytzg/2  (40) 

Repeating  this  process  for  the  particle  originatingat  hc  +  oz 
leads  to, 


Y2  =  Syt(z0+oz)/2  (41  ) 

Clearly,  when  Y2-Yi  =0y  ,  equation  (40)  subtracted  from  equation 
(41)  produces  another  expression  for  the  growth  of  lateral 
dispersion  defined  as, 

o  „ 


S„to../2 


(km) 


(42) 


IX.  Results  and  Conclusions 


Results 

The  simplest  and  quickest  method  of  comparing  the  model's 
output  with  the  expressions  presented  by  WSEG-10,  Norment,  and 
Bridgman  is  to  plot  all  four  simultaneously  for  a  particular 
value  of  shear.  Shown  in  Figure  20  is  such  a  plot  where  the 
shear  value  is  0.2  km/hr-km. 

The  WSEG-10,  Bridgman,  and  Norment  expressions  are 
plotted  exactly  as  they  are  shown  in  Figures  16,  17,  and  18.  As 
a  reference  point,  the  value  of  og  at  24  hours  is  included  on 
the  right  edge.  Again,  these  three  expressions  plot  as  straight 
lines. 

The  results  of  the  model  are  also  plotted.  The  shear 
contribution  is  extracted  from  the  data  developed  by  the  model 
for  the  total  lateral  spread  (as  shown  in  Table  XII)  in  the  same 
manner  that  was  applied  to  WSEG-10's  equation  (equations  1  to 
3). 

On  this  scale,  it  appears  that  the  model  predicts  a  os 
which  asymptotically  approaches  the  plot  of  os  predicted  by 
Bridgman's  equation.  Although  the  plot  of  WSEG-10's  expression 
intersects  the  model's  plot  between  3  and  6  hours,  it  obviously 
diverges  rapidly  from  the  model's  prediction  at  later  times.  The 
predicted  growth  from  the  WSEG  equation  is  sufficiently 
different  from  the  results  predicted  by  the  model  that  the  WSEG 
equation  will  be  disregarded  from  this  point  on.  Both  Norment 


and  Bridgman's  equations  will  be  examined  on  a  larger  scale 
before  discarding  one  or  the  other. 

Figure  21  shows  the  results  from  the  Bridgman  and  Norment 
equations  and  the  model’s  results.  Shear  remains  at  0.2  km/hr- 
km.  Additionally,  both  0yU  and  Oy^  (from  Table  XII)  are  averaged 
to  form  a  single  oy  for  the  model.  Averaging  the  two  oy  from  the 
model  causes  the  asymptotic  plot  of  Figure  20  to  become 
essentially  a  straight  line  parallel  to  the  Bridgman  plot  for 
this  time  domain.  Norment's  equation  produces  a  line  that 
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shear  values  are  examined.  All  other  parameters  are  the  same  as 
for  the  0.2  km/km-hr  case  just  examined.  Figure  22  displays  the 
calculated  values  for  both  Bridgman's  and  the  model's 
predictions  for  shear  values  of  0.6  and  1.0  km/km-hr. 

Conclusions 

This  study  has  examined  the  processes  that  cause  fallout 
particles  to  be  spread  in  a  crosswind  direction.  An  improved 
model  of  the  initial  stabilized  cloud  and  a  discrete  altitude 
dependent  fall  model  have  allowed  these  processes  to  be 
accurately  modeled.  From  the  widths  of  the  crosswind  patterns  as 
determined  by  the  model,  it  is  clear  that  none  of  the  initial 
three  equations  accurately  predict  the  crosswind  spread.  Dr. 
Bridgman's  equation  is  the  best  analytical  method  for  predicting 
the  kinds  of  og  that  the  model  found.  As  is  evident  in  Figures 
21  and  22,  it  is  valid  for  the  range  of  shear  values  currently 
considered  to  exist  in  the  lower  atmosphere. 

The  model's  disadvantages  are  the  computation  time  and 
the  unevenness  of  the  product  due  to  the  discrete  mesh  of  the 
particle  distribution  with  altitude.  The  disadvantage  of  Brid¬ 
gman's  equation  is  that  an  average  fall  time  must  be  assumed. 
This  assumption  leads  to  predictions  of  os  that  are  more  accu¬ 
rate  than  the  other  equations  analyzed  but  still  not  as  accurate 
as  the  model.  However,  the  close  agreement  between  the  new 
equation  formulated  by  Bridgman  and  the  results  of  the  model, 
make  an  excellent  case  for  using  Bridgman's  equation  for  predic¬ 
ting  the  lateral  growth  of  fallout. 


Second,  a  variety  of  particle-activity  distributions 
could  be  used  to  describe  different  initial  stabilized  clouds. 
These  could  then  be  matched  with  actual  cloud  data  with  the  goal 
of  determining  the  most  appropriate  particle  distribution  to  use 
in  conjunction  with  the  model  developed  here. 

Third,  an  analysis  of  the  <vz>  term  in  Bridgman's 
derivation  would  allow  more  detailed  study  of  the  validity  of 
his  assumption.  Furthermore,  it  could  lead  to  an  expression  for 
lateral  growth  that  is  a  function  of  weapon  yield. 

Lastly,  a  study  of  the  mu lti-burstscenario  could  be 
attempted.  Whereas  the  single  burst  can  be  modeled  with  a 


gaussian  distribution,  the  multiple  burst  cloud  and  fallout 
pattern  would  require  a  different  distribution. 


DELFIC  100  Equal-Activity  Group  Stabilization  Data 


Group 

Radius  (pm) 

hc  (km) 

c2  (km) 

1 

1917 

0.0 

0.206488 

2 

1064 

0.0 

0.987598 

3 

782 

0.0 

1.24583 

4 

629 

0.836141 

1.38594 

5 

529 

2.87594 

1.47751 

6 

457 

4.34459 

1.54344 

7 

403 

5.44608 

1.59289 

8 

361 

6.3028 

1.63135 

9 

326 

7.01673 

1.6634 

10 

297 

7.60827 

1.68996 

1 1 

272 

8.11822 

1.71285 

12 

251 

8.54658 

1.73208 

13 

232 

8.93414 

1.74948 

14 

216 

9.2605 

1.76413 

15 

202 

9.54608 

1 .77695 

16 

189 

9.81125 

1.78886 

17 

177 

10.056 

1.79984 

18 

167 

10.26 

1.809 

19 

157 

10.464 

1.81816 

20 

148 

10.6476 

1.8264 

21 

140 

10.8108 

1.83373 

22 

133 

10.9535 

1.84014 

23 

126 

11.0963 

1.84655 

24 

119 

11.2391 

1.85296 

25 

113 

11.3615 

1.85845 

26 

108 

1  1.4635 

1.86303 

27 

103 

11.5655 

1.86761 

28 

98.4 

1  1.6593 

1.87182 

29 

93.8 

1  1.7531 

1.87603 

30 

89.5 

1  1.8408 

1.87997 

31 

85.5 

1 1.9224 

1.88363 

32 

81.7 

12.0 

1.8871  1 

33 

78.1 

12.0734 

1.89041 

34 

74.7 

12.1427 

1.89352 

35 

71.4 

12.2101 

1.89654 

36 

68.4 

12.2712 

1.89929 

37 

65.5 

12.3304 

1.90195 

38 

62.7 

12.3875 

1.90451 

39 

60.1 

12.4405 

1.90689 

40 

57.5 

12.4936 

1.90927 

Group 

Radius  (pm) 

hc  (k®) 

0  z  (km) 

41 

55.1 

12.5425 

1.91147 

42 

52.9 

12.5874 

1.91349 

43 

50.7 

12.6323 

1.9155 

44 

48.6 

12.6751 

1.91742 

45 

46.6 

12.7159 

1.91925 

46 

44.7 

12.7547 

1.92099 

47 

42.8 

12.7934 

1.92273 

48 

41.1 

12.8281 

1.92429 

49 

39.4 

12.8628 

1.92585 

50 

37.7 

12.8975 

1.9274 

51 

36.2 

12.9281 

1.92878 

52 

34.7 

12.9587 

1.93015 

53 

33.2 

12.9893 

1.93152 

54 

31.8 

13.0178 

1.93281 

55 

30.5 

13.0443 

1.934 

56 

29.2 

13.0708 

1.93519 

57 

28 

13.0953 

1.93629 

58 

26.8 

13.1198 

1.93739 

59 

25.6 

13.1443 

1.93848 

60 

24.5 

13.1667 

1.93949 

61 

23.5 

13.1871 

1.94041 

62 

22.4 

13.2096 

1.94141 

63 

21.4 

13.23 

1.94233 

64 

20.5 

13.2483 

1.94315 

65 

19.5 

13.2687 

1.94407 

66 

18.7 

13.285 

1.9448 

67 

17.8 

13.3034 

1.94563 

68 

17 

13.3197 

1.94636 

69 

16.1 

13.3381 

1 .94718 

70 

15.4 

13.3523 

1.94782 

71 

14.6 

13.3687 

1.94856 

72 

13.9 

13.3829 

1.9492 

73 

13.2 

13.3972 

1.94984 

74 

12.5 

13.4115 

1.95048 

75 

11.8 

1 3 . 4258 

1.95112 

76 

11.2 

13.438 

1.95167 

77 

10.6 

13.4502 

1.95222 

78 

10 

13.4625 

1.95277 

79 

9.45 

13.4737 

1.95327 

80 

8.89 

13.4851 

1.95379 

82 


*.  -*.  -r. 


Radius  (ym) 

hc  (km) 

°z  (km) 

8.35 

13.4961 

1.95428 

7.82 

13.507 

1.95477 

7.32 

13.5172 

1.95522 

6.83 

13.5271 

1.95567 

6.35 

13.5369 

1.95611 

5.89 

13.5463 

1.95653 

5.45 

13.5553 

1.95694 

5.02 

13.5641 

1.95733 

4.6 

13.5726 

1.95771 

4.19 

13.581 

1.95809 

3.8 

13.589 

1.95845 

3.41 

13.5969 

1.9588 

3.04 

13.6045 

1.95914 

2.68 

13.6118 

1.95947 

2.32 

13.6191 

1.9598 

1.97 

13.6263 

1.96012 

1.62 

13.6334 

1.96044 

1.27 

13.6406 

1.96076 

.904 

13.648 

1.9611 

.473 

13.6568 

1.96149 

Appendix  B 


Data  for  Activity  Distribution  Plots 


3  hours: 


Radius 

Grounded  Location 

Acti vity-Agi 

( ym) 

(km) 

(Ci) 

85.5 

3.45673 

4.51 982E-4 

81 .7 

3.76633 

9.59174E-4 

78.1 

4.0525 

1 .44708E-3 

74.7 

4.31519 

1 .68805E-3 

71.4 

4.56274 

1.621 46E-3 

68.4 

4.78215 

1.35534E-3 

65.5 

4.99054 

1.01 696E-3 

62.7 

5.18892 

6.991 37E-4 

60.1 

5.35708 

4.58880E-4 

57.5 

5.53205 

2. 8674  IE-4 

12  hours; 


34.7 

12.8071 

3.66634E-4 

33.2 

14.3793 

9.02291E-4 

31.8 

15.7713 

1 .44682E-3 

30.5 

17.0255 

1.69530E-3 

29.2 

18.2353 

1 . 56075E-3 

28.0 

19.307 

1  .  19573E-3 

26.8 

20.3353 

7.80603E-4 

25.6 

21.3185 

4.4651 0E-4 

24.5 

22.1979 

2.45446E-4 

24  hours; 


23.5 

25.3690 

4.00393E-4 

22.4 

29.0203 

1 .03766E-3 

21.4 

32.2468 

1.57847E-3 

20.5 

34.9277 

1 . 68178E-3 

19.5 

37.7679 

1 .32815E-3 

18.7 

39.9412 

8.99450E-4 

17.8 

42.2687 

4.91 404E-4 

17.0 

44.2221 

2.53800E-4 
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